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Abstract 

Light scattering has proved to be one of the most powerful techniques for probing the properties of particulate systems. The 
purpose of this paper is to review the status of elastic light scattering. The emphasis is on recent developments, rather than being 
over-repetitive of earlier reviews, but considerable background is included with the aim of making the paper self-contained. 
There is an extensive summary of theoretical treatments, including both historical work and new. On the experimental side, 
while there is some discussion of well-tried methods, the emphasis is on recent techniques for measuring other properties as 
well as size. This includes, among others, the fractal treatment of agglomerates, determination of particle shape and measure- 
ment of refractive index. The discussion is broad rather than deep to provide a wide-ranging review of an extremely active field. 
© 1998 Elsevier Science Ltd. All rights reserved. 
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Nomenclature 

Roman tetters 



A 

c 

C 
D 
D f 
E 

A 

F(B % <t>) 
h 
H 
H 

m 

I 
/ 
I 

lm 

h 

J, 

k 



L 
m 



expansion coefficient in series for held 
cross-sectional area (m*) 
dimension less vector amplitude 
expansion coefficient in series for field 
speed of light (m s~'): mass concentration 
(kg m~ y ) 

cross-section (nr) 
diameter (m) 

diffusion coefficient (nrs^ 1 ); fractal dimension 
vector electric held (Vm* 1 ) 
volume fraction 

dimension less scattering function 

scalar Greens function 

height of ray above axis (m) 

spherical Hankel function of first kind of order n 

affinity exponent 

cylindrical Hankel function of first kind of order n 
vector magnetic held (A m" 1 ) 

irradience (W m~ 2 ) 
identity tensor 
intensity (W nT 2 ster" 1 ) 
imaginary part 

spherical Bessel function of order n 
cylindrical Bessel function of order n 
wavenumber (2ir/A ) (m~ 1 ); Bohzmann's constant 
(JK' 1 ) 

coefficient (m" 1 ); Kubelka-Munk absorption 
coefficient (m _l ) 

premultiplier in fractal expression for agglomer- 
ate 

cut-off of Gegenbauer spectrum 
length fm> 
refractive index 



Mi 
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a 

N 

p 

<? 

Q 

r 

f 

R 

R u R 2 



5 
t 

T 

m 

u 

V 

V 
W 

W 0 
x 



element of Stokes matrix 

second moment of size distribution 

number of particles per unit volume in elemental 

size interval dD (m~ 4 ) 

unit vector in the direction of n 

number of particles per unit volume (m~ ? ); 

number of particles in fractal agglomerate 

normalized phase function 

power (W) 

scattered power per unit volume (W m ~ 3 ) 
scalar flux (W m~ 2 ): 2k sin(0/2) <nT') 
efficiency 

radius of individual sphere in agglomeration 
radial position 

radial position at internal point 
radius (0/2) (m) 
Fresnel reflection coefficients 
radius of gyration <m) 

element of amplitude scattering matrix; Kubelka- 
Munk scattering coefficient (m _1 ) 
Poynting vector or radiative flux (W m~ 2 ) 
time (s) 

temperature (K); transmissivity 
Chebyshev polynomial 
T matrix 

particle velocity (m s~') 
gas velocity (ms* 1 ) 
volume (m*); visibility 
width (m) 
beam width (m) 

particle size parameter (ttZVX) 



Greek letters 



eCr) 
V 



permittivity of free space (Fm* 1 ) 
relative permittivity of panicle 
cos<e) 
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0 scattering angle 

£ M angular position of nth rainbow . 

A wavelength (m) 

Af fringe spacing (m) 

p viscosity (N sm" 2 ) 

permeability of free space (Hm" 1 ) 

v frequency (s" 1 ) 

£ characteristic cluster size 

p density (kg m" 3 ) 
pv, p H depolarization ratios 

t turbidity (KL) 

<f> azimuthal scattering angle; phase 

0 scalar field 

^ shape factor 

o> 0 albedo 

B solid angle (ster) 

Subscripts 

0 original or incident 

abs absorbed 

ae aerodynamic 

diff diffracted 

eff effective 

ext extinction 

geom geometrical optics 

inc incident 

pr radiation pressure 

sea scattered 

t transmitted 

ve volume equivalent 

Superscript 

* complex conjugate 

1. Introduction 

In numerous cases there is the need to understand the 
nature of mult i -phase systems. These occur in combustion, 
for example, whether particles are burning, being formed or 
are added as tracers. However, liquid drops and solid parti- 
cles have practical implications in very many industrial and 
everyday situations. 

In order to determine the properties of particulate systems 
and understand their behaviour it is necessary to acquire 
diagnostic tools. Paramount among these is electromagnetic 
radiation since it can be used remotely and, at the intensities 
normally used, does not interfere with the system under test. 
The advantages of visible wavelengths are not only that the 
radiation can be seen, but the wavelength is small 
(» 0.5 u.m) yielding excellent resolving power. 

When light interacts with a phase boundary a number of 
processes can occur, including reflection, refraction and 
diffraction. When the boundary encloses a discrete volume 
or when it has structure within itself, these processes are 



3 

often lumped together under the title "scattering". Analysis 
of the scattered light reveals the properties of the structure. 

The study of light scattering by small particles has a long 
and honourable history. The books by van de Hulst [ 1 } and 
Kerker (2] are classical reviews of the subject. More recent 
general works include the paper by Jones [3] and the books 
by Bohren and Huffman [4] and Bayvel and Jones (5]. 
Laser-based techniques for industrial applications have 
recently been discussed by Black et al. [6] Instrumentation, 
with particular emphasis on combustion applications, was 
reviewed by Jones [7], 

The purpose of this paper is to consider recent develop- 
ments in the theoretical treatment of scattering and its use 
for particle characterization. The intention is to be general 
without being too repetitive of the earlier reviews. Thus, 
instrumentation is not discussed unless it has some particu- 
lar novelty. However, while the aim is to emphasize recent 
work, an attempt is made to make the paper self-contained 
and it is hoped that sufficient background material is 
provided to make it understandable. 

The review begins with a definition of the basic variables 
used in describing light scattering, and then moves on to 
general theoretical treatments. These are of two types: 
those which may be considered as "rigorous" in the sense 
that, in principle, they can be calculated to any desired 
accuracy and those which are approximate. The former 
includes both analytical and numerical techniques which 
are limited only by computer power. The latter covers a 
wide range of theories, each of which is valid only over a 
limited range of parameters but which provide easy means 
of calculation and much insight. The review of theoretical 
methods concludes with a discussion of multiple scattering 
which occurs in dense media. 

Experimental methods for the measurement of particle 
properties are described in the second half of the paper. 
These cover techniques applied both to clouds of particles 
and to individual particles to measure parameters as diverse 
as size, concentration, refractive index, shape, structure and 
velocity. 



2. Bask definitions of scattering variables 

2.7. individual panicles 

When light is incident on a particle two processes can 
occur, it can be deflected or absorbed. For particles that 
are large compared with the wavelength of the radiation, 
the deflection processes are commonly described by reflec- 
tion, refraction and diffraction, though these will not 
describe certain phenomena such as the glory which arises 
from waves in the surface of the particle. For particles that 
are of the order of the wavelength in size or smaller, we 
cannot distinguish these processes and can only refer to 
"scattering". However, in general we group all the 
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Lens 



<b) 



Fig. I. Diffraction by an object: (a) Frcsnel diffraction; (b) Fraun- 
hofer diffraction. 



deflection phenomena and absorption under the heading of 

scattering. 

Diffraction and scattering arc often used synonymously 
and it is worth spending a little time to show how they differ 
and to what extent they are the same. We recognize that 
diffraction is usually divided into two classes, as indicated 
in Fig. 1. The first, as in Rg. 1(a), is FresneJ or "near-field" 
diffraction. In this scattered light interferes at the observa- 
tion plane with thai unscattered. Clearly, the scattered and 
diffracted light patterns are not the same. 

In the second, the observation screen is moved to a 
distance far from the scatterer — ideally to infinity. This is 
Fraunhofer or "far-field" diffraction as illustrated in Fig. 
Kb). Since parallel lines meet at infinity, the situation can be 
simulated by observation in the focal plane of a lens. Here 
the incident light is focused to a spot at F, whereas the light 
scattered through some angle 0 is focused at some other 
point such as O. The scattered and incident light waves are 
quite separate and do not interfere with each other. It follows 
that the diffracted and scattered light patterns are identical. 

Thus in the case of Fraunhofer diffraction the difference is 
only one of terminology. However, "diffraction" is usually 
used for an approximation valid for large particles in which 
a scalar wave interacts with their cross-section. The term 
"scattering" is reserved for vector wave interaction with 
the body of a particle. 

We now consider a single particle illuminated by radia- 
tion with irradience (incident intensity) / 0 (W m" 2 ). In all 
that follows we shall assume that the incident radiation has 
the form of an infinite plane wave. Practically this implies 
that the amplitude and phase must be sensibly constant 
across the particle. Certain of the results, such as the optical 
theorem (see below), can become invalid for shaped beams 
when the particle size is a significant fraction of the beam- 
width. This is an important aspect of the Gaussian output 
from lasers, for example. 



The total scattered power P^, is proportional to /,>, the 
constant having the dimensions of area, i.e. 

where C»* is the scattering cross -section. Similarly, for 
absorption 

Cat* being the absorption cross-section. Both of these 
processes remove power from the incident radiation, and 
their combined effect is extinction. The power extinguished 

is 

and 

C c xt -~ Cjcj + Cjrt* 

is the extinction cross-section. 
The power incident on the particle is 

where A is the cross-sectional area. The ratio of scattered to 
incident power is 



inc 



A w 



where Q ra is the scattering efficiency. It follows that 

The light scattered by a particle diverges over the surface 
of a sphere and its propagation must satisfy the wave equa- 
tion in spherical coordinates. In free space far from the 
particle the scattered electric field has the form of an 
outgoing spherical wave with an angular modulation. We 
write 



Since the field must remain transverse to the direction of 
propagation, which in this case is r, can only have terms 
in A t and A+. Thus 

From Maxwell's equations it can be shown that the 
magnetic field is 

in 



• - -(*) (M - U) 



/*0 

The Poynting vector, which determines the energy flux in 
the wave, is 

5= iRe(£x//*) 

which becomes 
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Fig. 2. Polarization states relative to the scattering plane. 

where F<0.<£) = |A(0.^)| 2 is the dimensionless scattering 
function. Since rr = I, the scattered power is given by 
integration over a sphere at radius r as 



sea 



= * r r 

*gJo Jo 



F<0,<f>)sin0d0d<£ 



and the scattering cross-section is 



C 



*f>Jo Jo 



F<0.d>)sin 0d0d* 



We note that 



^ = -J^<0,*) 

is the scattered intensity (W m~ 2 ster"'). 

If the scattering is isotropic, in the sense that it is inde- 
pendent of angle, then 

Using this we define the normalized scattering function, or 
phase function, by 



Pi0,4>) 



so that 

r2w 



4ttJo Jo 



l 



The extinction cross-section is obtained from the optical 
theorem [I, 2] 

When measuring scattered light, it is usual to allow the 
detector to rotate around the scattering centre in a plane as 
shown in Fig. 2. It is convenient to describe the incident and 
scattered light such that the polarization is resolved into 
components parallel and perpendicular to this plane. The 
fields are related through 

/%ca\ e^7*2 

which defines the amplitude scattering matrix. We note that 



for an isotropic particle, such as a sphere, no cross-polariza- 
tion is introduced and S 3 = S 4 = 0. Also, for any particle 
there can be no indication of the azimuth angle <f> at 0 = 0°. 
Then it also follows that, at 0 = 0*. 5 3 = S 4 = 0 and that 5, = 
S 2 = 5(0). Thus 



c Uc ° r 



and A = S(0)l Finally. 

C M1 = ^ lm[S<0)] 

There is one other cross-section which should be 
mentioned in passing, namely that for radiation pressure, 
Cp,. The rate of momentum of a light wave is given by 
(P/c& For light scattered through the angle 0 the rate of 
momentum in the original direction of propagation is 
changed to {PuJcfi. Thus there is a force on the particle 
given by 



sot 



(I -cos0)£ 



where cos 0 is the average of cos 0 obtained by integration 
of the scattered intensity over all values of 0. In addition, 
any power the particle absorbs stops the radiation and 
contributes a force given by (P^tcfa The total force is 



Thus the force per unit irradiance is 

T = ~[C„, - coT&C^Jz 
'o c 

By analogy with the previous definitions we define Cp, as 



Cpf — Cgji cos 0C|£j 



Since pressure is force per unit area, the radiation pressure is 
given by fip/c. 

Another property of the amplitude scattering is that if the 
incident and scattering directions are interchanged then 



which is the reciprocity theorem. This important relation- 
ship has been derived, for example, by Saxon (81. 

When we examine scattered intensity we need a means of 
describing all the polarization states and their relationships. 
This is achieved through the Stokes' matrix which relates 
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four % 'intensities* \ defined by 

U = £|£l + £ A .£j 
through 









l 








V 



II M|2 Af |3 Af |4 \ 
A^2i A/in ^24 
i*#3, Afi 3 Af M M M 



The Afy form the elements of the Stokes* matrix. They are 
given in terms of the 5, of the amplitude scattering matrix by 
Bohren and Huffman (41, for example. 

Then? are various symmetry relationships between the 
elements In particular it can be shown that, for a 
cloud of randomly aligned particles (1,9], 

(Mm 



Af I2 



0 
0 



0 
0 



0 
0 



0 \ 
0 

Mu I 

For isotropic particles, M u = and A/33 — ^44* If either of 
these inequalities is unsatisfied anisotropy is present. 

2.2. Particle clouds 

Particles very rarely, if ever, exist singly. They exist in 
clouds containing very large numbers. The actual number of 
interest to light scattering depends upon the circumstances, 
such as the volume of the test space as seen by the detector. 
Thus at low concentrations and with a small volume we may 
be looking at a single particle, but the presence of the others 
cannot be ignored in all cases. For larger volumes and 
concentrations we may be looking at many millions of parti- 
cles. As we shall see later, measurement techniques natu- 
rally fall into these two classes: single-panicle counting or 
cloud methods. 

The extension of single scattering theory to clouds can be 
performed at two levels, the simplest being at low concen- 
tration. The complications arise at high concentrations 
where multiple scattering and particle interaction effects 
become significant. 

For very low concentration the simple definitions of the 
previous section can be extended to a cloud in a straightfor- 
ward manner provided three conditions are met: 

). The panicles are randomly positioned and are of suf- 
ficient number that incoherent superposition may be 
applied. The random positioning means that the phases 



Table 1 

Maximum concentration of particulates to avoid direct interaction 
effects 



Particle radius (11m) 


Concentration (nT*) 


0.1 


to* 


1.0 


10* 


10.0 


10" 


100.0 


10»° 


1000.0 


10 7 



of the scattered waves arriving at the detector are also 
random. If the number is sufficiently large that the phases 
sum to zero, then intensities can be added directly. 

2. There is negligible multiple scattering. This means that a 
photon once scattered has a very high probability of 
leaving the cloud without encountering another particle. 
We shall examine the requirements for this shortly, but, 
approximately, it means that the transmission by the 
cloud must be greater than 60%. 

3. There is no interactive scattering. This implies that the 
electric fields should not interact directly. If this happens 
then there is a coherent structure to the multiple scatter- 
ing and the particles tend to behave to some extent as an 
assembly. In the limit they are touching and they behave 
as a single panicle or agglomerate. To avoid interactive 
scattering the evidence is that a separation of three 
diameters is thought to be sufficient [2]. The implications 
in terms of panicle size and concentration are shown in 
Table 1. The maximum volume fraction that can be 
occupied by panicles of any size is approximately 3%. 

If all three conditions apply we simply add intensities. 
Thus, if there are W* particles per unit volume all having 
the same size £>. the total scattered power per unit volume 
of space is 

P scslN — NP sea 

and, on dividing by the irradience. 



***** _ N p < 



k sca 



where X*. is the scattering coefficient. The absorption may 




Fig. 3. Transmission through a slab. 
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be similarly defined, and the extinction coefficient is 

More usually, there will be a distribution of particle sizes. 
We allow there to be n{D) 6D particles per unit volume in 
the sire range D to D + dD. Then 

iV = n(D) dD 



and 



' VC3 



CUnfD) d£> 



and so on. 

We now consider transmission through a slab containing 
particles, as in Fig. 3. The volume of a section of length dx is 
A dt. so that the power lost is 

-dP CM = IK^A dx 

or the loss of irradiant flux is 

-d/ = IK ext dx 

On integration 

/, = /„ exp(- £ K txt dx) 

If the cloud is uniform throughout the length, then 
/, = /„ exp(- K^L) 

The product K t JL is called the turbidity, r. 

The optical mean free path (the average distance a photon 
will travel between scattering events) is \IK W . This suggests 
that to avoid multiple scattering requires 

/,// 0 = e"' =0.368 

The actual value will depend upon the particle type. For 
example, if a particle is very heavily absorbing, so that 
scattering is relatively small, multiple scattering may not 
occur until interactive effects become appreciable. 

To have some idea of the limits to concentration we may 
examine large spheres. As we shall see later, for large pani- 
cles Q txy 38 2 so that 



irD 2 



Thus, for example, 10 ujti particles will require NL < 1 .6 x 
10 ,0 m~ 2 . 100 jun particles require NL < 1.6 x 10 s m" 2 . 
These concentrations are quite modest, and multiple scatter- 
ing is very common. 

A brief discussion of the treatment of multiple scattering 
is given in Section 4. 

3. General scattering theory 

We now briefly summarize means for calculating the 
expected scattering from particles. Fundamentally, there 



are two approaches. In the first, expansions are written 
down for the incident, internal and scattered waves and a 
solution is found for the expansion coefficients by applying 
the electromagnetic boundary conditions. In the second, the 
solution is found from an integral over the volume of the 
scatterer. In both cases there are limitations and complica- 
tions. 

J. /. Boundary condition method 

To obtain an analytical solution by this method, the fields 
are expanded in a coordinate system in which the tangential 
components of the field are readily obtainable. Thus for a 
sphere, for example, spherical coordinates are chosen. The 
tangential fields then lie in the 6 and <t> directions. 

In order to obtain a form for the expansion the wave 
equation must be separable in the chosen coordinate system. 
Further, there must exist orthogonality relations to enable 
calculation of the expansion coefficients. These require- 
ments limit the number of possible solutions to six coordi- 
nate systems: rectangular Cartesian; circular, elliptical and 
parabolic cylindrical; spherical; and conical |10 r 11]. 

The infinite cylinder has been studied in some detail as it 
is the most straightforward geometry. The solution for the 
case of normal incidence was first obtained by Lord 
Rayleigh [12]. The more general case of oblique incidence 
was solved by Wait [131, and is discussed by Kerker [2] and 
Bobren and Huffman [4J. Other solutions for radially stra- 
tified cylinders and cylinders with inhomogeneous shells are 
reviewed by Kerker [2], By considering the limits as the step 
in thickness and refractive index to go to zero, Kai and 
d'Alessio [14] were able to consider very large numbers 
of shells with high computational efficiency. Their program 
could easily accommodate up to 80 000 shells. 

A suggested scheme for the treatment of a finite cylinder 
was given by van de Hulst [I]. In this, the result for the 
infinite case was multiplied by the diffraction pattern of a 
slit of length equal to that of the cylinder. More recently, this 
approach was revisited by Wang and van de Hulst [15]. 
They claim good agreement in comparison with microwave 
measurements on models. 

The more useful, and most widely applied, solution is that 
for a sphere. This can be applied to real particles, though it 
must be emphasized that most particles are non-spherical. 

The solution for the sphere was first derived by Lorenz 
[16] and later independently by Mie [17]. It is discussed in 
detail by van de Hulsl [ 1 ], Kerker [2] and Bohren and Huff- 
man [4], among others. The elements of the amplitude scat- 
tering matrix are given by 

£\ n{n + 1) 



*»-l|nT>k* + *««l 
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Fig. 4. Scattering diagrams and efficiency factors: (a) intensity versus angle for* = 0.01 and m = 1.33: (b) intensity versus angle forjr *= 5 and 
m = 1 .33; (c) intensity versus angle for x = 20 and m - 1 .33; (d) intensity versus angle for x = 20 and perpendicular polarization, showing the 
effect of absorption; (e) extinction efficiency showing the effect of absorption; (f) efficiency factors for absorbing particles: (g) extinction 
efficiency in detail, illustrating the presence of singularities form = 1 .33. 
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The coefficients a* and b> for a non-magnetic homogeneous 
panicle have the form 



a n = 



= ztii!\z) are 



fo(i«r)fc'(*) - mfc 
K fc(ifu)t'W " **. '<«*)£,<*) 

where x = <Mz) = Vr,U) and f*(z) 
Riccati-Bessel functions. 
Also 

tt„(cos0)=« ^ P J< C0S 



The scattering and extinction efficiency factors are 

G iM = 4if2«+lMI«.Mb.f] 



06 



a „ = -jRc5;(2 n + i)K + 6»] 



There are several algorithms for the calculation of a„ and 
b n . These are discussed by, for example, Bayvel and Jones 
(5) and Bohren and Huffman [4], A FORTRAN program is 
listed in the back of Bohren and Huffman's book. It is found 
that the series are convergent and can be terminated after a 
certain number of terms. M Bohren and Huffman use 

It is evident that as the particles become larger more terms 
are required and. eventually, the computer time becomes 
excessive. Ultimately it becomes more efficient to use 
geometrical optics. 

Some typical results for scattered intensity and the 
efficiency factors are shown in Fig. 4. There are a number 
of typical points. For a very small particle with* 1 (Rg- 
4(a)) the scattering pattern is that of a simple dipole. For 
polarization perpendicular to the plane of measurement the 
scattered intensity is independent of angle and the scattering 
is isotropic. For polarization parallel to this plane there is a 
cos 2 6 variation. As the particle becomes larger, forward 
scattering increases more rapidly than that at other angles 
and quickly becomes dominant. The shape of this lobe for 
large particles (s4u,m at visible wavelengths) is well 
described by Fraunhofer diffraction theory. This shows 
that the height of the forward lobe is proportional to D A 
and its width to 1/D. The scattering pattern also becomes 
convoluted with a complicated lobe structure. The number 
of lobes is of the order x\m\. Absorption by the particle 
damps the lobe structure away from the forward direction. 
Again, for large particles where Fraunhofer diffraction 
holds, calculation shows that the forward scattering is not 



sensitive to refractive index. This is very important in prac- 
tice to particle sizing if refractive index is not known. 

The efficiencies initially rise rapidly with particle size, 
but then peak and go asymptotically to a constant with a 
decaying oscillation. Lock and Yang [18) have demon- 
strated that the oscillations are caused by interference 
between near forward transmission and diffraction. 

For the extinction efficiency the value of the constant is 
two. This would suggest that the particle will block off twice 
the light falling upon it, an effect called the "extinction 
anomaly". The explanation lies in the fad that two 
processes are occurring: diffraction and the geometrical 
optics effects of reflection, refraction and absorption. The 
efficiency for each of these processes is one, so that 

For very large particles the diffraction pattern is very 
narrow. All real detectors, including the eye, collect light 
over a finite angular range. Thus the diffracted light is 
collected and is not measured as a loss. The measurement 
yields 



as would be expected. In order to ensure that the true extinc- 
tion efficiency is measured, it is necessary to restrict the 
angular collection range of the detector to be much less 
than the width of the forward lobe. Thus a measure of trans- 
missivity is not necessarily the same as a measure of extinc- 
tion, and great care is required. 

If the efficiencies are calculated in fine detail, with very 
small steps in x, a number of extremely sharp peaks and 
troughs are observed. These are scattering resonances asso- 
ciated with the zeros of the denominator in the expressions 
for a n and b„. The expansion coefficients become infinite at 
these points. 

The extension to radially stratified spheres is quite 
straightforward and is discussed, for example, by Kerker 
[2], A FORTRAN program to calculate scattering from 
two concentric spheres is listed by Bohren and Huffman 
|4], Multiple concentric spheres have been considered by, 
for example, Johnson [19] and the application to radially 
inhomogeneous spheres has been discussed by a number 
of authors (e.g. Perelman 1201). 

Kai and Massoli [21] considered the limits as the steps in 
size and refractive index go to zero and obtained a simple 
and speedy algorithm. This method has enabled a very large 
number of shells to be considered and various aspects of 
radial inhomogeneity have been investigated. An example 
is the influence of a radial variation of refractive index on 
* the resonances in Q txt . Kai and d' Alessio [22, 23] and Kai el 
al. [24] concluded that their location depends primarily on 
the average value of the retractive index near to the surface. 
If the refractive index only changes in a thin layer near the 
surface, then an equivalent homogeneous sphere can be 
found. 
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Fig. 5. Primary rainbow angles for 10 urn droplets during hearing 
Un = 1.43 - 0.05/1) and cooling <m = r J3 + 0.05*). FSSM refers 
to a parabolic profile of refractive index using the finely stratified 
sphere model. The Mie theory results are for homogeneous spheres. 
t After Kai ei aJ. (251-1 

Another area of interest lo potential diagnostics is the 
position of the rainbow. Kai et al. 125] used the stratified 
sphere model to indicate the ways in which this position 
fluctuates. Figure 5 compares the variation of the rainbow 
angle from a homogeneous sphere with one having a vari- 
able profile of refractive index. In this case a homogeneous 
drop with a refractive index of 1.43 enters a hot environ- 
ment. As the drop heats up a temperature profile develops. 
Initially only the outer edge is heated but then the higher 
temperature penetrates further until, eventually, the drop is 
once again homogeneous at the final temperature where the 
refractive index is 1.33. As the drop cools the process is 
reversed. In the stratified sphere model the refractive 
index is allowed to vary over a given shell of increasing 
thickness with a parabolic profile. In the Mie calculations 
the sphere is homogeneous with a linear variation of refrac- 
tive index. Jt is evident that there is considerable difference 
between the two models and that significant error could be 
introduced by the assumption of homogeneity. 

The detection of inhomogeneity via measurements of the 
scattering matrix was investigated by Bhanti et al. [26] They 
used a stratified sphere model to predict elements of the 
matrix for coatings of soot on both water drops and coal 
particles. 

The case of a sphere containing an eccentric inclusion has 
been considered by Borghese et al. [27, 28] and Videen et al. 
(29] In a similar approach. Fuller |30. 31] considered the 
coated sphere in terms of multiple scattering with a view to 
developing a model for spheres with multiple inclusions. 
Not only does the presence of an inclusion affect the absorp- 
tion cross-section very considerably, but it is strongly 
dependent upon the position and orientation relative to the 
incident wave. This is illustrated in Fig. 6. 

Videen et al. f32] and Chylek et al. [33) also looked at the 



case of carbon particles enclosed within water drops because 
of the relevance to absorption of radiation in the atmos- 
phere. The former authors suggested the use of effective 
medium predictions and concluded that the extended meth- 
ods for the calculation of equivalent refractive index 
removed any restriction on the size of the internal grains. 

The Lore nz- Mie theory can also be used to calculate the 
internal fields of spheres. This has become of interest parti- 
cularly because of studies in laser heating, and also because 
of interests in radiation pressure. Examples of such studies 
are the papers of Lock and Hovenac [34] and Lai et al. [35] 

Notable among solutions for other shapes is that for the 
spheroid due to Asano and Yarnarnoto [36]. This theory is 
complicated and computationally expensive. To date scat- 
tering has only been calculated for quite small spheroids. A 
recent advance due to Kurtz and Salib [37] has modified the 
formulation of the boundary conditions to yield a more 
convergent form. Another approach due to Martin [38] is 
to apply a perturbation to Mie theory for the case of slightly 
deformed spheres. Coated spheroids have been considered 
by Farafonov et al. [39] and Voshchinnikov [40]. 

As mentioned earlier, the theory outlined above is for 
infinite plane waves. However, with the advent of the 
laser, attention has been drawn to scattering out of Gaussian 
beams. Early work in this area was due to Tarn and Corri- 
veau [41]. This problem has been vigorously pursued by 
Gouesbet and co-workers in recent years. Although the 
general problem is very complex* they have been able to 
develop a quite simple form for sufficiently large panicles 
on the beam axis by employing van dc Hulst's "localization 
principle" |1, 42-45). A rigorous justification of the loca- 
lization principle for both on-axis and off-axis beams has 
been propounded by Lock and Gouesbet (46, 47]. Ren et al. 
1 48) have discussed die evaluation of the beam shape 




a 



Fig. 6. Specific absorption cross-sections for a 0.1 |un carbon 
sphere eccentrically located in a water drop of diameter 5 tun. 
The refractive indices are 1.5 - H>.5 and 1.33, respectively. a x 
and a> arc the radii of the water drop and the inclusion, and d t* 
the radial distance of the inclusion. The cross-sections have been 
averaged over polarization. (After Fuller (3t|.) 
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coefficients for laser sheets. Lock [49] and Ren et ai. [50] 
have proposed improved algorithms for scattering by a 
sphere out of a Gaussian beam, and Gouesbet [51] has 
examined scattering by an arbitrarily located and oriented 
infinite cylinder. Hie validity of this generalized Mic theory 
has been demonstrated by Guilloteau et al. [52] by observing 
light scattering by particles levitated on laser beams. 

Interesting new developments occur when this theory is 
applied to situations where the beam width is of the order or 
less than the diameter of the particles. For example, Lock 
(53] demonstrated mat a narrow Gaussian beam incident 
near the edge of a particle can produce rainbows up to 
ninth order. Also, Hodges et al. [541 find that diffraction 
theory fails if the beamwidth is less than the diameter. 
Beam confinement also affects extinction and this has 
been explored by Gouesbet et al. [55], who found that addi- 
tional terms were needed in the optical theorem. Lock 156] 
has shown that when the beamwidth approaches infinity, 
Q nx -* 2 in agreement with current knowledge. However, 
when the beamwidth approaches zero. &*t — 1 because of 
the loss of the diffraction term. 

Attention is now being given to beam shapes of a more 
general nature [57-59]. Lock and Hodges [60] have inves- 
tigated scattering of an off-axis non-Gaussian beam by a 
sphere, and Onofri et al. [61] have performed similar calcu- 
lations for a mulrj -layered sphere. 

Most particles in nature are not spherical. Indeed most 
particles have shapes which cannot be strictly treated by the 
boundary condition method. When faced with this situation 
we have to be pragmatic and accept that an analytical solu- 
tion is not possible. There are two routes we can lake. We 
can either approximate or we can use a numerical method. 
The latter is not an approximation in so far as the calculation 
can, in principle, be continued to any desired accuracy. The 
only limitations are computer time and space and numerical 
precision. Nonetheless, certain assumptions are often made 
in setting up the numerical scheme which can make it only 
approximate, or limit its application within certain bounds. 

There is a set of numerical schemes based on the bound- 
ary condition method. The principle here is to presume that 
the fields can be expanded in terms of some basis functions 
which enable them to be written down and matched at the 
boundary. Perhaps the simplest is the point matching 
method in which it is assumed that the series can be termi- 
nated after N terms. The consequence is that there are then 
W unknown coefficients. Selection of a sufficient number of 
points on the surface enables a set of simultaneous equations 
to be written down which can be solved for the unknowns. 

The question is that, if it is this straightforward, why is it 
not more widely used? The problem lies with the assump- 
tion that the series is convergent. Of course, it is if the 
matching takes place over the surface of a simple shape 
for which a series solution is known, but it is not necessarily 
so for other shapes. There has been considerable discussion 
on this over the years (e.g. Lewin [62]). 
Joo and Iskander [63] modified the point matching 
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method by expanding the field using spherical harmonics 
at various points along the axis of an extended object. In 
each zone they used point matching but, additionally, 
ensured that the fields in different regions agreed where 
they overlapped. They also used a least-squares fit for the 
coefficients and claimed coasiderable reductions in compu- 
ter time. 

A similar concept is the multipole technique (e.g. Ludwig 
(64]) which uses a set of point sources within the particle, 
the total field of which approximates the true field. A linear 
combination of such points is found which satisfies the 
boundary conditions, and matching at a discrete number of 
points reduces the problem to a matrix inversion. The matrix 
can be over-determined and the method of least squares 
applied. Generally, this avoids the problem of ill-condition- 
ing of the matrix. A further extension of this idea is mode 
matching (65]. A generalized multipole technique has been 
developed by Al-Rizzo and Tranquilla [66] which enabled 
them to examine quite large elongated and composite 
objects. 

An attempt to avoid the convergence difficulty lies with 
the extended boundary condition method (EBCM). Water- 
man [67] considered a wave incident on a perfect conductor. 
This generates electric currents in the surface which cancel 
the internal field. By the use of analytical continuation it was 
shown that the field cancellation held throughout the pani- 
cle's interior. The scattered field is obtained from an integral 
of these currents over the particle surface. The method has 
been extended to dielectric particles by, for example, Barber 
and Yeh [68]. 

The number of terms, N, of the series expansion required 
for convergence increases with size and refractive index, but 
also with deviation from spherical shape. For any particle it 
is possible to draw inscribed and exscribed spheres of radii r\ 
and r e . The spherical harmonic expansion can be shown to 
be convergent at radii r < r, and r > r e . The Rayleigh 
hypothesis is that the series is also convergent within r, < 
r < r e . It is this assumption which is dubious. The EBCM, 
which was supposed to be exact, makes a similar assumption 
and suffers from the same kind of defects. Lewin (62] 
showed that this is due to hypersensitivity to minute errors 
of the field on the matching circle. If the true shape differs 
significantly from a sphere then difficulties will arise with 
convergence. This can make the method very time-consum- 
ing for all but the smallest, only slightly perturbed sphere. 

Numerous calculations have been made for small parti- 
cles where the shape is essentially spherical with a small 
perturbation. Most popularly, the shape is described by an 
expansion in Chebyshev polynomials such that 

r = R[l +5T B (cos0>] 

and 6 is small. For example, Wiscombe and Mugnai [69] 
compared Chebyshev particles and spheres. They found 
reasonable agreement between the two types for panicle 
size parameter, jc, less than three to five. After that the 
differences grow rapidly. A universal feature was that the 
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Fig. 7. Light scattering polar diagrams of irregular particles. The 
solid tine is for an ideal sphere of the same size, (a) A rough and 
slightly elongated glass sphere (50 pin), (b) A fly ash particle 
( 15 |tm, m * 1.55 + iO.OI ). (After Killiager and Zerull [70].) 

lobe structure in the scattering pattern became averaged out 
for the non-spheres. They also found that panicles with 
concavities scattered substantially more than totally convex 
ones. Differences in intensity could be summarized as 
below: 
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An example is shown in Fig. 7(a). This was measured 
experimentally by Killinger and Zerull (70) for a randomly 
oriented, slightly irregular, 50 ujn diameter glass sphere. 
Very small deviations from spherical shape cause distinct 
differences in the scattering properties. The conclusions are 
largely in agreement with those of Wiscombe and Mugnai. 
They are not universal, however, as witnessed by Fig. 7(b) 
which Killinger and Zerull obtained for an absorbing fly ash 
particle of 15 ujn diameter. 

Later. Mugnai and Wiscombe [71 ] attempted to represent 
scattering from irregular particles by that from a size distri- 
bution of spheres, which has the similar effect of averaging 
out the lobe pattern. However, the results were not easy to 
interpret universally. The authors commented that even after 
many hours on a Cray supercomputer they did not have 
sufficient results to come to any firm conclusions. 



There have been several attempts to combine the EBCM 
with other methods to improve speed and convergence. For 
example, Iskander et al. [72] proposed an iterative proce- 
dure starting with a known approximate solution. They 
claim that this approach is more stable and quicker, and 
used it to study elongated particles, such as spheroids [73, 
74]. For this, the fields were expressed as expansions in 
spherical harmonics about several points along the axis of 
the particle. The EBCM has also been applied to ellipsoids 
by Schneider and Peden [75] and to finite cylinders by 
Ruppin [76] and Kuik et al. [77) The latter authors comment 
that for random orientation the solution does not converge to 
that of an infinite cylinder as the axial ratio increases, though 
it does for single aligned cylinders. They also state that 
ellipsoids do not model finite cylinders well because of 
the sharp edges. Barton [78] has studied layered Chebyshev 
particles illuminated with focused Gaussian beams. 

The EBCM is usually combined with the r-matrix 
method and the two terms are used interchangeably. The 
r-matrix simply relates the expansion coefficients of the 
scattered field to those of the incident field, i.e. 



where EBCM is used to compute the elements of [71 This 
matrix contains all the information about the scatterer, such 
as size, shape and refractive index. Rotation of the scatterer 
can be allowed for simply by rotating the incident field, 
which changes a wJix and b VM . A discussion of this method 
and of its use for the determination of particle shape has 
been give, for example, by Barber and Hill [79]. 

Because of the intensive computational requirements 
coupled with the influence of rounding errors on the ill- 
conditioned matrix which has to be inverted, there has 
been an upper limit to the particle size parameter that can 
be calculated, of the order 25. Mishchenko and Travis [80] 
have managed to extend this by up to 2.7 times by the use of 
extended precision computing. 

Mishchenko (81. 82] has considered the problem of light 
scattering by randomly oriented particles with axial sym- 
metry. A method is described to generate the components 
of the Stokes' matrix directly for assemblies of such parti- 
cles which avoids the need to perform many individual 
calculations. Storage requirements are reduced and it is 
claimed that computation times can be reduced by a factor 
of up to 20. Contour plots of elements of the scattering 
matrix for polydispersions of randomly oriented spheroids 
have been provided by Mischchenko and Travis [83]. Orien- 
tatkmal averaging in the r-matrix method has also been 
discussed by Khlebtsov [84]. In a further paper by 
Mishchenko [85], simplifications are discussed for the 
calculation of Q nt for ellipses. 

A general review of the application of the r-matrix 
method to non-spherical panicles has been provided by 
Mishchenko et al. [86]. 
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An extension of numerical techniques to non-spherical 
particles illuminated by a Gaussian beam is discussed by 
Barton and Alexander 187]. They provide a numerical 
scheme for solving the boundary conditions, and the results 
are expressed in terms of expansions in spherical harmonics. 
The method was tested for off-axis spheres and ellipsoids. 
Khakd et al. [88, 89] have treated the same problem by 
expanding the Gaussian beam as a spectrum of plane waves. 

Recent work has concerned the extension of the r-matrix 
method to interactive and multiple scattering (90-93). 

3.2. Integra! equation methods 

The alternative numerical methods are based on the inte- 
gral formulation of the scattering problem. It has been 
shown by various authors that a solution of the scalar 
wave equation is given by 

*<r> = + *5 f yi W) - \W)G(r, ?) dV 

where r and r' are two points in space. G(r, ?') is the appro- 
priate Greens function. In infinite cylindrical coordinates 

In spherical coordinates 

c *.|?-?'| 



47r|r - r'| 



^ toc (r) is the incident field at? The boundary conditions are 
automatically taken into account and the radiation condition 
is satisfied. 

If r is a point internal to the particle, the integral equation 
can be used to calculate the internal field. If f is an external 
point, then it can be used to calculate the scattered field 

*Un = $lM*')-W')GCr.7)4V' 

External to the particle, it is usual for practical purposes that 
|r| » |r'| and the far-field approximation to the Greens func- 
tion may be used. For spherical coordinates this is 



■\k^r cos 0 



4ur 



G{rJ') 
where 

cos 0 m cos 0cos 0* + sin 6 sin & cos(<£ - <t>) 

It can be seen that there are two difficulties with these 
equations. The first is that the unknown field appears inside 
the integral. The second is that there is a singularity in the 
Greens function when |r - r r \ = 0. 

Despite these difficulties the integral formulation has 
formed the basis of a number of calculation methods for 
the scattered field. Some of these are approximations 
where a form is assumed for the internal field. This is the 
subject of the next section, but the most commonly used 



assumption is that 

= 

which is the Bom approximation. Once the internal field is 
given, the scattered field can be calculated without the 
problem of the singularity arising. 

Another possibility is to expand the internal field by a 
known set of functions (commonly spherical harmonics) 
with unknown expansion coefficients. Then an attempt is 
made to match the coefficients, perhaps by some method 
such as least squares. 

The most common numerical method is to subdivide the 
internal volume into a number of much smaller volumes. An 
assumption is then made about the nature of the field within 
each sub-volume. Typically, if the sub-volume is suf- 
ficiently small, the internal field is taken to be an unknown 
constant. Then we have 

The problem then reduces to solving N linear simultaneous 
equations. There is still, however, the problem of the singu- 
larity. This can be tackled in two ways. One is to draw an 
exclusion zone around the point and leave this small volume 
out of the problem. This zone is then made successively 
smaller and the solution is found in the limit as the volume 
of the zone tends to zero. This is equivalent to taking the 
principal part of the integral. 

The other way out is to let the volume around the point be 
a sphere. Since we assume that the field is constant, we may 
restrict the calculation such that r, is at the centre of the 
small volume. It can be shown that on integration the singu- 
larity disappears. The assumption of spherical shape for this 
one small volume should contribute very small error. The 
accuracy can be improved if the singularity is removed by 
using an expansion in Taylor's series to two terms [94]. 

For the purposes of electromagnetic scattering the vector 
equation is needed. This is (95] 

x ( i+ ~k* 'v'j<Xrf)W 

where I is the identity tensor. An alternative form is given 
by Saxon 196J which in some ways is more tractable. This is 

£(?) = + *6 \ yt W) - l]£(r')C<?.r') dV 

As we are dealing with three-dimensional space, any vector 
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can have three components. As a consequence, when the 
method of division into N sub-volumes is used, we end up 
needing to solve 3N simultaneous equations. 

Perhaps the simplest approach to discretization of the 
integral is to make the sub-volumes simple cubes or rectan- 
gles (e.g. Richmond [971; Hagmarm et al. [98]). This is the 
basis of the method of moments on which there is a consid- 
erable literature. A good discussion of this general area is 
given by Newman and Kingsley [99]. 

An improvement that has been proposed is to use terra- 
hedra rather than cubes because these are belter able to fit a 
curved surface [100, 101]. Hagmann et al. [98], Schaubcrt ct 
al. [100) and Mendes and Arras [102] also propose 
improved sets of basts functions, rather than simple 
constants, which permit the use of larger sub-volumes 
and, thus, reduce computer demand. It has been shown by 
Lu and Chew [103] that computational efficiency can be 
improved by noting that a group of small volumes can be 
represented by a set of surface sources on their common 
boundary. This fact is used to continuously replace the inter- 
nal sources with equivalent larger volumes and results in a 
slight decrease in computer time. 

A logical extension is to allow all the sub-volumes to be 
spheres. Indeed in the limit as the radius of each sphere 
approaches zero, they behave as simple dipoles. This 
forms the basis of the coupled dipole approximation in 
which the dipoles are considered as equivalent to the mol- 
ecules making up the particle. 

The coupled dipole method probably began with the 
paper by Purcell and Pfennypacker [K04]. It has been 
explored particularly by Singham and Bohren [105, 106], 
who developed a method of successive approximations 
rather than solve a large number of simultaneous equations, 
which they call the "scattering order method". The range of 
validity has been extended by Mulholland et al. [107] who 
included the first magnetic dipole. This enabled the primary 
sphere diameter to be increased from 0.06 u.m to 0.12 jxm. 
Similarly, Okamoto (108] expressed the dipole polarizabil- 
ity of the primary spheres as the first term (a,) of the Mie 
series expansion, and achieved calculations for size para- 
meters up to 1 .5, or about 0.24 \xm in the visible. An algo- 
rithm for improving the speed of the calculations has been 
given by Venizelos et al. [109] Draine and Flateau [1 10] 
have reviewed the method. 

The iterative approach requires an equivalent medium 
approximation which allows for the interactive Meld 
between the dipoles and alters their polarizability, or refrac- 
tive index. A number of authors have recently been explor- 
ing this. Ku [111] has compared various models and 
concluded that the one provided by Iskander et al. (1 12] is 
the most satisfactory. Wolff et al. (1 13] prefer (he extended 
Bruggeman theory proposed by Rouleau and Martin 
[ 1 14, 1 15]. Yet another proposal from Draine and Goodman 
(1 16] allows for the radiation reaction of the dipoles and 
dispersion in the medium. 

The most appropriate model may depend upon the 



volume fraction of the dipole particles. Lumme and Rahola 
1117] give 

"in ~ 1 ~ m" - 1 

+ 2 + v(ml n - 1) /v m> + 2 + - 1) 

For v = 0 this reduces to the Maxwell -Gamett equation and 
for v — 1 it is equivalent to Bruggeman theory. From numer- 
ical calculation these authors conclude that 

/v » 

0.52 -0.3 to 0.3 

0.41 0 
0.3 0.3 

Lakhtakia and Vikram (1 18| have also proposed a modi- 
fled form of the Maxwell -Gameu equation which allows for 
the finite size of the particles. This takes the form 




where 



0 _ (m^wiO 2 -! 

1 -[(«2/«i) 2 -l]{|(I -iton f )e*~« - l] 

Subscript 1 refers to the surrounding medium and subscript 
2 to the particles. The authors claim that the method is 
suitable for \kam } \ £ ir/5» where j is 1 or 2. and 0 < /„. < 0.2. 

Drolen and Tien [1 19] and Dobbins and Megaridis [120] 
have explored treating spherical agglomerates approxi- 
mately as solid spheres with an equivalent refractive index 
calculated using Maxwell -Gamett theory. Dittman [121] 
has also proposed a method for agglomerates employing 
an equivalent refractive index but only including first- 
order multiple scattering, which is claimed to be sufficient. 

The successive approximation method is efficient for 
small particles but does not converge for large particles. 
As an example, it is necessary that x < 3 for m = 1.33. x 
increases as m decreases. For the full coupled dipole 
approach very many dipoles may be needed. For example, 
Ckxxlrnan et al. [ t22J found for m m 1 .33 + iO.Ol and x = 10 
that 3 X 10 4 were required. The number scales as (jr|m|>\ 
The computational requirements quickly become excessive, 
even using a Cray-2 [123]. In their review Draine and 
Flateau [1 10] show that the method can be accurate for 
\m\ < 2 and they discuss a solution technique which permits 
up to I0 5 dipoles and jt s 10. 

A simplifying suggestion has been proposed by Kumar 
and Tien [124] which allows the internal fields of all the 
particles to be the same except for a phase term due to their 
relative positions. This approximation removes the need to 
solve simultaneous equations and is claimed to have good 
accuracy. 

Methods for improving computer efficiency in solving the 
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Fig. 8. Ratio of efficiency and volume equivalent size parameter, 
x rtI for a cluster of 40 primary spheres of size parameter x r m = 
1.6 + i0.6. (After Mackowski 1133].) 



simultaneous equations have also been discussed. For exam- 
ple, Flaieau et a!. [125] showed that for a rectangular array 
of sub-volumes the matrix has a block-Toeplitz structure. 
They were able to exploit this to improve the efficiency of 
the computation. Yung [1261 suggested certain reductions 
from consideration of the requirement of reciprocity. Inter- 
estingly, however, Charalampopoulos and Parugrati [1271 
have raised the possibility that reciprocity may not hold 
for agglomerates, perhaps as a consequence of multiple scat- 
tering. 

Bourrely el al. [1281 nave suggested that computational 
efficiency can be improved by not using a uniform distribu- 
tion of dipoles. In their approach the particle was divided 
into inner and outer regions. The inner region contains 
larger dipoles with greater separation. Small closely spaced 
dipoles are only used near the surface. They claim their 
method to be more accurate, quicker and space-saving, 
and that it is better at describing the particle shape. 

Recent applications of the coupled dipole method include 
porous particles [113, 129, 130], finite cylinders [131] and 
particles on a surface [132]. 

The accuracy of the method has been explored by Mack- 
owski [133] by comparison with exact theory for multiple 
sphere clusters. Examples for the extinction and absorption 
efficiencies are illustrated in Fig. 8. The coupled dipole 
approach tends to underestimate the increase in absorption 
by a cluster. In a later paper Mackowski [134] showed that 
the internal field can be non-uniform, leading to errors in the 
coupled dipole method. However, an electrostatic analysis 
was explored for the case where the size of the cluster is 
much less than the wavelength, which was able to resolve 
the problem. Also, Rouleau [135] compared the calculations 
for dense clusters of up to 30 spheres against the exact 
solution. It was concluded that the dipole model was not 
satisfactory for large particles, and that equivalent refractive 
index models (such as Maxwell -Garnett) may be prefer- 
able. 

Michel [ 1 36] and Michel ct al. [ 1 37] have explored scat- 



tering by individual irregularly shaped particles using the 
so-called "strong permittivity fluctuation theory'*. They 
begin with an equivalent medium description and then 
extend it using a perturbation method. Essentially they 
solve the wave equation in the form 

V X V X E(f) - J^ ff <r)£(r) = k 2 [e{f) - e tff (r)]E(r) 

where e(r) - s eff (r) is small. The statistical properties were 
obtained by a Monte Carlo method to generate the clusters. 
For a cluster of 64 grains each of size parameter jc = 3.3 X 
10" 5 . they claim a scattering enhancement of twice that 
from an effective medium theory alone. 

The coupled dipole method and the method of moments 
have been reviewed and compared by Lakhtakia and 
Mulholland [138]. They demonstrate that the two methods 
are equivalent. In a similar vein, Flateau et al. [139] have 
compared the two techniques for a system of two spheres 
and have found good agreement for x *» 10. 

A method gaining in popularity is the forward difference 
time domain (FDTD) technique. In this, the electric and 
magnetic fields are obtained by stepwise integration of the 
wave equation across the particle. Schlager and Schneider 
[140] have produced a wide-ranging review and Yang and 
Lion [141] give a good description of the method 

With advances in computing, new techniques promise 
radical improvements in speed and handling ability. In 
this context, Cwiket al. [142] have discussed the application 
of massively parallel computation to integral equation tech- 
niques in light scattering. 

Reviews of numerical techniques have been given, for 
example, by Bates 1143], Yaghjian [144], Miller [145] and 
Bohren and Singham [146]. 



4* Approximations 

Rigorous solutions of the scattering problem, whether 
analytical or numerical, are complicated and tend to involve 
considerable computer time and storage for all but the 
simplest shapes and smallest particles. For certain ranges 
of size and refractive index, approximations may be used. 
Not only are these simpler and more readily computable, but 
they often provide insight into the scattering process. 

We shall consider the various approximations roughly in 
terms of increasing size. 

4.1. Rayleigh scattering fx < /; x|m - /| <£ /; 

The limitations imposed ensure that the incident electric 
field is uniform both within and around the particle. It then 
oscillates as a simple dipole. For incident polarization along 
the jc -direction it can be shown that 

= £o *o>(cos S cos 4>B - sin M>) 

where 4> is the angle between the plane of measurement and 
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x . A similar result is obtained for incident polarization along 
v, except that 4> is replaced by <f> + ir/2. For the case of 
orthogonal measurement planes 

(£K<7 :)(£) 

From this 

For unpolarized light 



Cam = -477*0 Refiy) 



The cross-sections are 

For a sphere 

= 3V m 2 - l 
? ~ 4ir/» 2 + 2 

where V is the volume of the particle. It follows that 



*r-l 



m 2 +2 



Oah, = 4jt Im 



V w2 + 2 / 



For ellipsoids there are three polanzabtlities correspond* 
ing to the three axes. These are given by 



V 

7j ~~ A- 



m 1 - I 



4tt 1 + (nr - \)Lj 

The Lj are related through 
L x + + ^ = I 

In general the are given by complicated integrals. 
However, simplification occurs for certain situations [1.2]. 
For a cloud of randomly oriented ellipsoids 

By comparison with Mie theory. Kerker et al. [147] have 
suggested that for agreement within 1%, Rayleigh theory 
can be used for x\m\ £ 0.2. For agreement to within 10%. 
x\m\ s 0.5. Similar work by Ku and Felske [148] suggested 
that Rayleigh theory was correct to within 1% for all jt sS I 
except for son panicles of low absorption (my < 1.2; m> < 
I). There seems to be a certain element of contradiction 
between these two conclusions. 

The theory can be extended to slightly larger particles by 
expanding the scattering coefficients in the form of polyno- 
mials in jr. Penndorf [149, 150] has provided an expansion 
for the efficiency factors of spheres, which is useful up to 
about x = 1.4 and w - 2 for dielectric spheres. For absorb- 
ing material the limit is x = 0.8 for 1.25 s m x s 1.75 and 
m 2 22 l. This expansion has also been discussed by Selamet 



and Arpaci [ 1 5 1 ]. A series expansion for the field scattered 
by ellipsoids has been given by Stevenson [152], which 
appears to be limited to kb s 0.95 where b is the semi* 
major axis. 

4.2. The Rayleigh-Gans-Debyc (RGD) or Born 
approximation ( jm - /| «: /; x|m 

The two titles of this approximation arise due to different 
approaches, but they produce the same result. The RGD 
concept can be understood with the aid of Fig. 9. Because 
the scatterer is "soft" (|m - 1| < I) there is very little 
reflection of the incident held. Also since (x\m — lj 1) 
there is very little phase shift inside the panicle. As a result 
of these two constraints, the internal field is approximately 
the same as the incident field in the absence of the particle. 

Since the incident field is divergenceless, the integral 
equation for the scattered wave becomes 

E^(r) = *6 J v , {m 2 - l)£ i0c (r')O(r,r') dv' 

For incident polarization along i it is found that 



and the intensities are 



(t)=^**(T) 



,//|.oCOS0 

Kerker [2] gives values for R(6, <f>) for a number of cases. 
Two examples are: the sphere 

= 3/, («y« 

and the infinite cylinder at normal incidence 
tf(0,tf» = 2/,(«yu 

where u — 2x sin(0/2) and j t and J% are spherical and cylin- 
drical Bessel functions. 

The range of validity for spheres has been investigated by 
Kerker et al. [153). Farone et al. [154] and Heller [155]. For 
accuracy to better than 10%, m s 1.25 and x — • 0 are 
required. The maximum value of m decreases as x increases. 




Fig. 9. II lustration of the RGD method. The origin of coordinates is 
at O and the scattered wave is received at P. The vndisturbed inci- 
dent wave is scattered by the element of volume dV. 
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Ft*, 10 Scattering by irregular particles in the RGD approximation: (a) influence of height of the irregularity for a mean size parameter °f *0* 
HI) compart of scaLng by irregular panicles (solid line) and a cloud erf spheres (dashed line) with the same mean size parameter of 10 
and a size distribution equal to the probability distribution of radius. (After Al-Chabbi and Jones [162].) 



Latimer and co-workers [156-158) have compared 
models for ellipsoids, and developed equivalent sphere rela- 
tions which can be used for m < 1 .5. Barber and Wang 1 1 59} 
compared ROD with EBCM for a prolate spheroid form = 
1.05 and 2x|m - 1| ^ 1 -0. They concluded that errors of up 
to 20% can occur depending on orientation. The worst case 
is when the major axis is parallel to the incident electric 
vector. Best agreement was found for random orientation. 

A loose agglomerate of small particles is a candidate for 
treatment by the RGD method, since it may be argued that 
the incident wave would propagate through the structure 
with relatively little change. This approximation removes 
the need to solve simultaneous equations and is claimed to 
have good accuracy. The range of validity of such an 
approximation for a fractal structure was explored by Farias 
et al. [160, 1611. Fractal structures are discussed later. 

A recurring problem in light scattering is the treatment of 
irregular particles. One method would be to calculate the 
scattering by a large number of particles individually and 
then add the results together allowing for random position 
and orientation. However, since the particles are random in 
nature, it is preferable to employ a statistical analysis on a 
cloud directly. This saves considerably on computer effort. 
Al-Chalabi and Jones [162] performed such an analysis 
using the RGD approximation. They found that the result 
depended upon the probability distribution of radius (mean 
radius and standard deviation) and a correlation function in 
the surface. As might be expected, increasing height of the 
irregularities caused increasing smoothing of the scattering 



polar diagrams and enhanced backscatter. There was also a 
suggestion that the irregular particles could be treated as a 
cloud of spheres with a size distribution equal to the prob- 
ability distribution of radius, provided that the height of the 
irregularity was not too large. This can be seen in Fig. 10. 

A drawback to the RGD approximation is that it does not 
permit the introduction of polarization effects due to aniso- 
tropic particles. Haracz et al. [163] used a modified version 
in an attempt to overcome this problem. In this, the basic 
Rayleigh-sized spherical sub-volumes were replaced by 
ellipsoids, which is equivalent to introducing an anisotropic 
refractive index. This approach was also used by Jones and 
Savaloni [164] to explore scattering by finite cylinders at 
random orientation. Ruppin 176) examined the validity of 
this calculation by using EBCM. Good agreement was 
found for axial ratios greater than 10, and. generally, for 
m £ 1.1 or for jr £ 1. 

4.3. The Wentzel-Kramers-Brillouin (WKB) 
approximation 

The RGD or Born approximation can be extended in a 
number of ways. A simple approach, for example, is to 
allow for the phase shift in the incident wave as it penetrates 
the particle. For a sphere we may then write 

£(r)-£oe^ (/f+ ^ 

where z f is distance along the axis. 
The Wentzel-Kramers-Brillouin (WKB) approximation 
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Fig. II. Illustration oftheWKB method. 

is similar except thai the phase is allowed to vary with 
distance from the axis according to the particle shape, 
presuming rectilinear propagation. Klett and Sutherland 
1165] have noted that this approximation is much improved 
if allowance is made /or interna) reflection from the back 
face of the particle. In this so-called °two-wave** WKB 
method, the internal field in a sphere for a ray entering a 
distance ft from the axis, as in Fig. 1 1 , would have the form 



where s - >JR 2 - h 2 . Jones et al. [166] have generated error 
contour charts for the two- wave method. They demonstrated 
that this was superior to both the RGD and single-wave 
WKB models. Nonetheless, the range of application is 
very limited especially at angles greater than 90°. 

Further extensions may include allowance for reflection 
loss on incidence by using the Fresnel equations, and for 
refraction. The eikonal approximation is similar [167-169]. 
Similar methods have been employed by Sharma and 
Somerford [170] to provide simple equations for backseat - 
tering by soft obstacles. Roy and Sharma [171] have 
compared various soft particle approximations. 

4.4. Fraunhofer diffraction fx > V; |m - /| » /) 

This well-known approximation is valid for large obsta- 
cles and unpolarized light and scattering close lo the forward 
direction. For a sphere. Hodkinson [172] gives 

The term in the square bracket is the Airy function. It 
suggests that the shape of the scattering pattern depends 
only on size and is independent of refractive index. Indeed, 
it results in the well-known expression 

sin*= — 

for the angular position of the first minimum in intensity. 

Jones f 173] has produced error contour charts comparing 
Fraunhofer diffraction with Mie theory. For transparent 
particles the discrepancy is less than 20% for m > 1.3 and 
.r a 20 3jun at visible wavelengths). For absorbing 
particles, however, errors of less than 20% occur for m t > 
1.2 and .t £ 6 I pun at visible wavelengths). 



For non-absorbing particles of very low refractive index 
(< 1.1) there is some error even for large particles. This is 
caused by transmission of the incident wave which inter- 
feres with that diffracted, leading to anomalous diffraction. 
Kusters et al. [ 174] concluded that serious errors could occur 
if this effect was not taken into account. However this is 
simple to do in practice. 

It is commonly found that forward scatter is also insensi- 
tive to particle shape (e.g. Pollack and Cuzzi [175] and 
Zerull et al. (176]). This was explored for irregular particles 
by Jones 1177] and Al-Chalabi and Jones [178], who 
confirmed the lack of sensitivity to shape. As discussed 
for the RGD approximation, the important parameters 
were the probability distribution of radius and the correla- 
tion function in the surface. The diffraction approach also 
suggested that the cloud of irregular particles could be 
modelled as a cloud of spheres with size distribution equal 
to the probability distribution of radius for small height of 
irregularity. A diffraction approach to scattering by irregular 
crystals has been discussed by Muinonen (179). 

Complicating factors are introduced by the nature of laser 
beams and the physical properties of optical components. 
Diffraction by a panicle in a Gaussian beam has been 
explored by Grehan et a). [180] and Lock and Hovenac 
[181]. Good agreement was found between measurement 
and calculation using the generalized Mie theory. As part 
of their study the latter authors demonstrated how the 
Kirchhoff integral can be derived from rigorous theory. In 
addition, lenses introduce phase changes. This effect has 
been pursued by Lebrun et al. [182], who applied thin lens 
theory to provide corrections to the diffraction pattern. 

If the concentration is too high, multiple scattering occurs 
and alters the diffraction pattern so that Fraunhofer theory 
no longer applies. However, empirical evidence suggests 
that, in practice, the extinction can be as high as 50% before 
multiple scattering has to be considered (183, 184]. 

Various attempts have been made to correct for multiple 
scattering. A correction at the onset of multiple scattering by 
use of a first order term from radiative transfer theory has 
been considered by SchnaWegger and Glatter [185]. 

More complete attempts have been made. The first of 
these was due to Felton et al. ( 1 86] with further development 
by Hamidi and Swithenbank 1 1 87] and Cao et al. [ 1 88] The 
scattering volume was divided into a series of slices, each 
sufficiently thin thai only single scattering occurs within it. 
The first slice sees only the incident wave and scatters some 
fraction F according to Fraunhofer diffraction. The second 
slice sees the incident wave reduced in intensity by F plus 
the light scattered by the first slice, and so on. By following 
this procedure through the whole volume, the light distribu- 
tion at the detector rings is found. The authors provided 
correction procedures which were claimed to extend the 
applicable range of the diffraction method to extinctions in 
excess of 90%. 

A different approach has been employed by Hirleman 
[189 1 which uses statistics to predict the small-angle light 
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distribution after many scattering interactions. His results 
agree well with those of the other authors except at high 
ring numbers (large angles). Hirleman attributes this to the 
fact that, in the previous method, a photon scattered out of 
the system of slices cannot be scattered back in, whereas his 

method allows for this. 

Woodall et al. [190] attempted to get away from plane 
media by examining the case of a cylinder, which is more 
suitable to the case of a spray. Their theoretical approach 
involved the use of the radiative transfer equation and 
Monte Carlo methods. They found that, due to multiple 
scattering, the errors in the parameters of a log-normal 
size distribution could be 25% for the Sauter mean diameter 
and 50% for the standard deviation for the case of t = 3. 

4.5. Anomalous diffraction (x » 1; |m - 7| <£ I) 

At very low refractive index, the particle transmits light 
almost without deflection. This then interferes with the 
diffracted light, producing anomalous diffraction. The scalar 
diffraction integral equation for a circular aperture becomes: 

where s' is the distance from the aperture to the screen or 
detector. K{&) is the inclination factor! 191] and 

8 = 2(m- \hj* 2 ~ ? 

For non-absorbing particles the extinction efficiency is 
m 4 sin a> 4(1 - cos (a) 

e- = 2--^- + — j — 

where <a = 2x(m - I). 

Chylek and Klett (192], Streekstra et al. \\91i\ and Four- 
nier and Evans 1194] have discussed non-spherical particles. 
In the latter case, the extinction efficiency for randomly 
oriented spheroids was checked for aspects ratios between 
0.5 and 4 with 1.01 £ m t ^ 2 and 0 ss /n 2 i I. More 
recently, Evans and Foumier [195] have provided an expres- 
sion for the direct calculation of the average extinction effi- 
ciency (&J for randomly oriented spheroids. Evans and 
Foumier [196] have provided means of directly calculating 
the averaged extinction efficiency factor. 

Sharrna [1971 suggests that a more appropriate criterion 
for validity of anomalous diffraction may be \m - l| 2 <£ 
\m + 1| J as opposed to \m - l| « 1. It was found by 
comparison with Mie theory that anomalous diffraction 
could work for m as large as 2. In later work Sharrna 
[198] has suggested an improved form for C„i obtained 
from expansion of the eikonal method. He proposes 

S(0)„ = C(m)S<0) ADA 

C(m) = 1 + i(m + l) 2 (m - 1X2 - mf 

This is claimed to be valid for xim 1 - l) 2 < 4. 



Stephens [ 1 99] performed similar calculations for circular 
cylinders. He found, surprisingly, that reduced to the 
correct form for the limit as x — 0 as well as for jt — oo. 

Chylek and Li [200] have discussed the boundary 
between the RGD and anomalous diffraction regions 
where x » I and x\m - l| 1. In particular, they derive 
expressions for the efficiency factors. Fot a complex refrac- 
tive index m = m, + vn 2 they obtain 

Q.» = VoniJ + * 2 [K " l) 2 -*!)i f J/ « 

where A is the cross-sectional area and V is the volume, I is 
the geometric path inside the particle, and 



-u, 



The extension of anomalous diffraction theory to allow 
for vector properties has been discussed by Bourrery el al. 
[201]. The method was applied to large spherical particles 
by an eikonal approximation to Mie theory. Anomalous 
diffraction by anisotropic scatterers has been considered, 
for example, by Meeten [202]. 

There have been a number of interesting recent studies. 
For example, Maslowska et al. [203) made comparison with 
the coupled dipole method for calculating scattering by 
cubes. They found that anomalous diffraction was not satis- 
factory for certain orientations and that care must be exer- 
cised in the application to non-spherical particles. 
Whitehead et al. [204] applied the technique to liquid crystal 
drops immersed in a polymer and Khkbtsov (205) to fractal 
clusters. 

4.6. Geometrical optics (x » /; x|m - /| > /,) 

Geometrical optics may be treated by a combination of 
Fraunhofer diffraction in the forward direction with reflec- 
tion and refraction at larger angles. Discussions and summa- 
ries are given by van de Hulst [1], Kerker [2], Emslie and 
Aronson [206] and Ungut et al. [207]. The latter authors 
compared the results with Mie theory for $ < 20° and 
found good overall agreement for x ^ 90 (s 15 \xm in the 
visible) and reasonable agreement for x £ 60 10 \ixn in 
the visible) at a refractive index of 1.5. However, agreement 
failed at certain specific sizes due to interference effects. 

Gtantschnig and Chen [208] produced an approximate 
geometrical solution in closed form. Agreement with Mie 
theory was good at small angles for D ^ 5 u.m in the visible. 
They also concluded that 3^ oc D l provided 

jrAot a 4ir 

where Ao was the angle subtended by the detector at the 
scattering centre. 
Shipley and Weinmann 1209] made comparisons with 
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Mic theory in (he range 200 < x ^ 4520 for m = 1 J33. 
They found very good agreement except at angles near graz- 
ing incidence and those associated with the rainbow or 
glory. 

There has been considerable discussion on the correct 
way to calculate efficiency factors. Cohen et al. [210] and 
Acquista et al. [21 1] propose 

Jim <?«. « I + £ # fl i? i| 2+ H 2 ] cos $s{n 6(16 

where Ri and R 2 are the Fresnel reflection coefficients. 
Senior [212] has been able to prove that 

lim Q ext = 2 

for perfectly conducting spheres. Nussenzveig and 
Wiscombe (213| have given 

- 2(. + 

Bohren and Nevitt (214] have produced an equation for 
£?«b* which appears to be accurate over a wide range of sizes 
and absorption indices. This is 

0 4 ™> 

(/n l + l) 2 -(m,-l) 2 e- 

where 

F= 1 + — J ml - 1, T = 4jdw^ 

■ 

It agrees very well in both the limits t < 1 and r — • oo. 

Kokhanovsky and Zege [215] have proposed equations 
for the efficiency factors due to spherical potydispersions. 
They essentially employ geometrical optics but include the 
influence of edge effects which are determined by subtract- 
ing the pure geometrical optics results from those of Mie 
theory. In a further paper [2161 they provide a general 
review of approximations for efficiency factors. 

There is continuing work on the application of geometri- 
cal optics to non-spherical particles. Ravey and Mazeron 
|217, 218). Hovenac [219]. Lock (220. 221] and Macke et 
al. [222] have produced theories for large spheroids. The 
latter authors suggest that ray tracing is more effective for 
spheroids than for spheres and can be used for x a 60 except 
for the prediction of polarization effects. Similarly, 
Mishchenko et al. [223] examined the application to finite 
cylinders by comparison with the 7*-matrix method and 
concluded that geometrical optics cannot adequately 
describe polarization effects. 

Bottlinger and Umhauer [224, 225] have considered the 
applications of ray tracing to individual irregular particles in 
which they were treated as a set of successive slices. Muino- 
nen et ai. (226] also looked at this problem but treated the 
panicles as random structures constructed from log-normal 



size distributions of height and allowing for the correlation 
function in the surface. 

Other interesting situations considered include inhomo- 
geneous particles (227], hexagonal crystals (228] and ice 
crystals with inclusions [229]. 



5. Multiple scattering 

5./. Independent scattering 

As we saw in Section 2, if the path length through a 
particle cloud is greater than or of the order of the optical 
mean free path, then multiple scattering becomes signifi- 
cant. This can be visualized by imagining streams of 
photons entering the scattering medium. If the particle 
concentration is sufficiently low, a photon once scattered 
has a low probability of being scattered again before leaving 
the medium provided the depth of the particle cloud is not 
too large. If the depth or concentration increases then so 
does the probability of further scattering events. In dense 
or extensive media a photon may be scattered a very large 
number of times before eventually emerging. 

As described above, the scattering events are independent 
of each other. It was suggested in Section 2 that a separation 
of three diameters may be sufficient for this. Provided this 
condition holds, the individual scattering events may be 
described by a theory appropriate to a single particle (e.g. 
Mie theory for a sphere). The only unknown is the amplitude 
of the incident wave, which must be obtained from a knowl- 
edge of the previous scattering events. 

Further, if there is a large number of particles randomly 
positioned in space, then incoherent superposition applies 
and intensities may be added directly. This is a stochastic 
process, and, indeed, one way of predicting multiple scatter- 
ing is the so-called Monte Carlo method [230-239]. As its 
name implies, this is a technique which uses random numbers 
to predict the direction of scattered photons, the probability 
of any direction being given by the appropriate scattering 
theory. The probability of distance travelled before a further 
scattering event is determined by the particle concentration, 
and a random number decides the actual distance. 

The Monte Carlo method is straightforward and simple to 
apply. The major drawback is in the computer time 
involved. According to Poisson statistics, the error after N 
scattering events is >/N, so thai the fractional error is \t\fN. 
This means that a very a large number of "tries** has to be 
made if the multiple scattering is to be predicted to any 
degree of accuracy. 

The Monte Carlo method employs successive singular 
scattering. There have been attempts to tackle this approach 
analytically. For example, Hottel and Sarortm [240] give 
expressions for the transmissivity and reflectivity of a 
plane slab. Similar concepts were employed by Deepak 
and Green [241) and Zahavi (242). At sufficiently high 
concentrations the random scattering process results in 
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Fig. 12. Illustration of terms in the radiative transfer equation. 

photons behaving in a manner similar to diffusion. This was 
the concept used by Rosseland [243]. 

Perhaps the most popular method for tackling multiple 
scattering is the radiative transfer equation with the multi- 
flux method of solution. This has the advantage of simplicity 
and the ease of physical understanding if the number of 
fluxes is small. The simplest is the two-flux method which 
leads to the Kubelka-Munk (244] equations which are 
widely used wherever there are reflecting layers. Typical 
examples are the paint and the food industries to describe 
thin films, such as of butter. The radiative transfer equation 
is also used particularly in radiative heat transfer, where an 
extra term is added to allow for thermal emission. 

To obtain the radiative transfer equation for multiple scat- 
tering the approach is similar to that in Section 2, except that 
we have to allow for radiation incident from other directions 
in the particle cloud being scattered into the direction of the 
incident beam. 

The situation is as in Fig. 12. Radiation of intensity 3' 
incident within the solid angle dft' is scattered through the 
angle ©' into solid angle dft in the forward direction. The 
incident intensity, as opposed to irradience, is 3. As 
previously, the scattering and absorption losses add to 



dP, 



3tf„ t dc dA dft 



Now, however, we have to add the gain due to in-scatter. For 
simplicity, we examine a single particle illuminated by the 
intensity 3'. The power incident is 



3'^dft' 



From the definition of the phase function, the fraction of this 
scattered through the angle ©' into dft is 

^-rte')Gsc dn 



For N such particles per unit volume the increase in power is 

dp 2 = y dfi'^p(e')f2« do n dv 

= 3Xe')<M'— dz<Wdft 
y ' An 

However, this is the contribution from one direction of inci- 
dence only. Integrating over all directions and adding the 
two terms we find the total change in intensity to be (since 
dP = d3 dA dft) 



f 3'p{e')dft' 



d^ — . t , ^sca 

— = -3K CJU + -P- 
dz 4ir 

We now make a modification to allow for the fact that we 
may be considering a ray not propagating normally to the 
faces of the slab. In general it may pass at an angle, as in Fig. 
13. Then the path length ds is 

ds = dz/cos 6 

and we have 

dZ 



cos 



Fig. 1 3. Illustration of a ray passing at the angle 6 through a slab. 



Writing 7f = cos 0, a>o = K^JK** where a> 0 is the albedo, 
and t = K nt x. where t is the turbidity, we find 

Kubelka and Munk [244] proposed the pair of equations 

^ - -{K + S)q+ + Sq- 
dx 

or 

for fluxes travelling in the positive and negative directions. 
From the radiative transfer equation it can be shown for 
isotropic scattering (p<6) — 1) that 5 83 and K = 
2/fafe. In general there is no such simple relationship. 
However, relationships have been found for a number of 
cases (e.g. Bayvel and Jones [51). 

In principle, the multi-flux method can be as accurate as 
required by taking many fluxes and expanding the phase 
function such that it resembles real particle scattering (e.g. 
Mudgert and Richards [245]). 

In view of the development of lasers there has been much 
interest in the transmission and reflection of collimated 
Gaussian beams and pulses. Examples of the former include 
the work of Beckett et al. [2461, Stotts [247], Look et al. 
[248], Zege and Kokhanovsky (249] and Ito [250]. Pulse 
spreading and broadening has been considered, for example, 
by Prahl et al. [251], van de Hulst and Kattawar [2521, 
Shifrin and Zolotov [253] and Brewster and Yamada [254]. 

Marquet et al. [255] used radiative transfer theory to 
predict the intensity profile of a transmitted beam, which 
they characterized in terms of the maximum and full 
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widih at half the maximum intensity. From this they were 
able to measure the scattering and absorption coefficients 
using a CCD array to display the profile. 

5.2. Dependent scattering 

When the particles approach closer than the limits 
suggested above there is direct interaction of the electric 
fields. The scattering by an individual particle is no longer 
independent since the incident field is modified by the 
presence of all the others. Also, the scattering becomes at 
least partly coherent and amplitudes have to be added. 

There are essentially two approaches to the problem of 
dependent scattering. The first is to develop the integral 
equation formulation in such a way that the volume integra- 
tion is carried out over the volumes of all the particles 
present. In its simplest form this is an extension of the 
coupled dipole approximation discussed in Section 3.2. 
where the dipoles are formed into groups each of which 
represents a particle (e.g. Katlawar and Humphreys [236]). 

The alternative approach is the boundary condition 
method. In principle, the method of point matching could 
be applied to solve the boundary conditions on every sphere. 
This would be very time-consuming and would suffer from 
the usual convergence problems associated with this 
method. 

To provide a rigorous solution it is necessary to be able to 
expand the boundary conditions on any one sphere about the 
centre of the coordinate system. This was achieved for two 
spheres, for example, by Levine and Olaofe [257] by the use 
of hispherical coordinates. More recent studies on two 
spheres include those by Barton et al. |258] and Fuller 
f 259 1, A study of a carbon particle both inside and attached 
externally to a water drop was studied by Fuller [260. 261] 
with relevance to absorption by atmospheric aerosol. 
Mishchenko et al. [91 1 have applied the /-matrix method 
to the problem of two spheres. 

More generally, a vector addition theorem has been 
developed which allows this procedure to be carried out 
for many spheres [262-270]. This approach has also been 
employed by Quimen and Kreibig [271] in the limit of very 
small particles, which then becomes equivalent to the 
coupled dipole approximation. Hamid et al. [272], Popovic 
and Notaros 1273] and Xu [274] have proposed iterative 
techniques which do not involve matrix inversion and. 
thus, avoid problems with ill-conditioning. 

The suggestion thai a separation of three diameters is 
sufficient to ensure that interaction between scattcrers is 
negligible has been explored by Ivezic and Menguc [275] 
for two Rayleigh-sized spheres. They found that a separa- 
tion of three radii was sufficient for interaction not to influ- 
ence absorption. However, for scattering the separation 
must be greater than twice the radius divided by the size 
parameter of the assembly. From their calculations they 
suggested that the ratio K<JK^ may be used to indicate 
the presence of agglomeration as well as the phase function. 



Other situations that have been examined include two 
infinite cylinders [276-278]. An experimental study of the 
case of two parallel fibres has been undertaken by Padrna- 
bandu and Bickel [279]. Arrays of cylinders have been 
explored by Bever and Allenbach [280]. Cooray and Ciric 
[28 1 ] give an analytical solution for two arbitrarily oriented 
dielectric spheroids and Sebak [282] discusses TM scatter- 
ing by two parallel elliptical cylinders. The problem of two 
spheres has been extended by Hamid et al. [283] to the 
situation where they are both stratified. Both the full bound- 
ary condition solution and an iterative technique were used. 



6. General comments on experimental methods 

.We have seen that the properties of the scattered light 
depend upon the particle size, shape, orientation, refractive 
index, and concentration. Thus, in principle, it is possible to 
use light scattering to measure alt of these parameters. This 
section reviews techniques aimed at measuring such proper- 
ties of particles. It will be assumed that visible wavelengths 
are employed, unless it is specifically stated otherwise. 

There are two approaches to light scattering measure- 
ments: 

1. integral methods. In this case large numbers of particles 
are viewed simultaneously. The scattered intensity is 
assumed to arise from incoherent superposition. This 
implies moderate to high concentrations and a significant 
sample volume. The lower the concentration, the larger 
the required volume. The resulting intensity comes from 
integration over space, and a spatial average is obtained. 
A spatial distribution can only be recovered directly if 
the total cloud is somewhat larger than the test volume. 

2. Particle counting. In this group particles are viewed one 
at a time. This implies moderate to low concentrations 
and a small sample volume. The higher the concentra- 
tion, the smaller the volume must be. The size distribu- 
tion is built up over time and a temporal average is 
obtained, which differs from the spatial average. This 
is easy to see. If we imagine two particle sizes with the 
same concentration but one having twice the velocity, 
then the faster particles will be counted at twice the rate 
and their concentration will be apparently doubled. If a 
spatial average is required then a conversion can only be 
made if the velocities of the panicles are known. If a 
spatial distribution is needed this must be built up over 
time by continuously repositioning the test volume. 

Particle counting introduces another problem. In the inte- 
gral techniques it is usually assumed that the particles are 
uniformly distributed within the test space. Any one small 
part of the sample volume sees the same size distribution 
and concentration as any other. It is then unimportant 
whether the light intensity varies within the test space. (Of 
course, to ensure that this is so, the test volume must be 
small in comparison to the scale of any variation within 
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the particle cloud.) In particle counting the position of the 
single particle is of paramount importance. A small particle 
near the centre of a Gaussian beam may yield the same 
scattered intensity as a large one near the edge. Thus there 
is a trajectory problem. 

There are three ways of overcoming this difficulty. The 
most obvious is to ensure that there is no variation of inten- 
sity within the test volume. This is the "top-hat profile" 
method. Unfortunately, the methods available for producing 
such a profile which also retain the coherent structure of the 
beam throw away much of the available light. The second 
technique employs deconvolution using previous knowl- 
edge of the structure of the test volume. Thirdly, it is 
often possible to choose a method of measurement involving 
relative variables (such as an intensity ratio or phase differ- 
ence) rather than the absolute scattered intensity. 

This highlights an important distinction between absolute 
and relative measurements. The absolute measurement of 
scattered intensity requires a knowledge of the incident 
intensity and of the response of the detection system. Thus 
the light source must be stable or monitored by a separate 
detector, and there must be no changes in intensity along 
either the path to the test volume or from there to the detec- 
tion system. The response must normally be obtained by 
calibration. A relative measurement has neither of these 
requirements in general, and is to be preferred if practicable. 
Typical relative measurements would be intensity ratio, for 
example at two or more scattering angles (including polar 
diagram methods or the shape of the forward lobe) or at 
various wavelengths, polarization ratio, and in laser Doppler 
methods the visibility and phase difference. Of these the 
visibility method can present a problem as it relies on an 
intensity balance between two beams. 

The advantage of using scattered intensity lies in its 
simplicity and its wide dynamic range, which is only limited 
to be between the noise level of the detection system and 
saturation. This can be extended by the use of neutral 
density filters and higher incident light intensity. Relative 
methods vary both in dynamic range and ease of use. The 
choice of appropriate technique depends upon the circum- 
stances, but the most widely used are reviewed later. 

Finally, it is worth noting some problems which 
commonly arise in the presence of gases, particularly hot 
gases. These include: 

1. Density gradients. The refractive index of a gas is in- 
versely proportional to its temperature and proportional 
to its pressure. Temperature and pressure gradients give 
rise to refractive index gradients which deflect the light 
beam (the "schlieren effect"). The linear displacement 
is proportional to the gradient and the path length 
through the medium. This may be considerable in 
large-scale systems. It can lead to uncertainty in the pos- 
ition of the test space, which may continually be in 
motion due to turbulence, spreading of the light, yielding 
uncertainty in incident and scattered intensities, and 



fluctuations in received intensity. A problem particular 
to anemometry is that the fringes in the test space may be 
continually moving, leading to errors in the measured 
velocity [283 J. 

2. Gases and particles. Gases are inevitable in most 
systems and particles are very common. Their signifi- 
cance to light scattering is twofold: they may interact 
with the incident and scattered light and they may radi- 
ate. The gases most commonly present (i.e. air) do not 
absorb in the visible and. apart from the effect of density 
gradients discussed above, they have no relevance unless 
particles of molecular size are being studied. However if 
infra-red or ultraviolet wavelengths are employed, for 
example, this is not so. Particles absorb, scatter and radi- 
ate at all wavelengths. The effects of absorption and 
scattering are to reduce the incident light intensity at 
the lest space and the transmission of scattered light to 
the exit window. With the path lengths typical of large- 
scale combustion rigs, for example, these effects can be 
very significant. They can be avoided by using relative 
measurements. Thermal radiation adds intensity at all 
wavelengths. Since the energy is spread over a broad 
spectrum, its influence can be greatly reduced by the 
use of narrow pass interference filters. Where appropri- 
ate, phase-sensitive detection can also be employed. 
Relative measurements do not necessarily help over* 
come the effects of radiation. 
3. Access. Many practical systems are enclosed and access 
is limited. This can dictate the technique to be employed. 
Often there are only two or three ports which severely 
restricts the possible scattering angles. If there is only 
one port then the use of backscattering is decreed. An 
alternative is to develop probes which can be inserted to 
lead the light in and out, but care must be taken to mini- 
mize disturbance to the system. Another problem with 
access is the presence of windows which inevitably 
become dirty. This reduces the incident and scattered 
light and, as before, presses the use of relative methods. 

The discussion below begins with integral methods and 
then moves to particle counting. In both cases the trend is to 
move from the smaller particles to the larger. Finally, there 
are sections devoted to the determination of properties other 
than size, such as refractive index and shape, to specialized 
techniques like tomography and inertia measurement and 
the new interest of particles on or near surfaces. 

7* Integral methods 

7. J. Rayleigh particles 

We begin by considering light scattering techniques for 
measuring the properties of particles which are very small 
compared with the wavelength, where the Rayleigh approx- 
imation holds. The intensity scattered by Rayleigh particles 
is too weak to detect from individuals, so that integral 
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scattering methods must be employed Perhaps the simplest 
method for absorbing panicles, such as soot, is to use ratio 
of scattering to extinction as proposed by van de Hulst [1]. 
We saw in Section 4. 1 that the scattering is proportional to 
volume squared while absorption is proportional to volume. 
Where the volume is very small this means that where 
absorption exists it dominates. Consequently 

Since 

oc AT?, oc N& 

where the bar denotes average over all sizes, it follows that 



'sea 



o 3 



Thus a size is measured which is biased by the size distribu- 
tion. The concentration can be recovered from the extinction 
measurement 
The ratio of scattering to extinction also involves 



m 

This shows that refractive index can be an important 
unknown. Its dispersion can also be significant. However, 
experiments by Bruce et aL [284] have suggested that oc 
I /A over a wide range of wavelengths. This variation was 
also found by Choi et at. [285], who measured extinction by 
the soot in a rich acetylene/air flame and obtained the 
geometrical properties using a gravimetric technique. 
They suggested 

and 

K t * 8.6 ± 1.5 

However, Koylu and Faeth [286] combined extinction 
measurement with an RGD analysis for aggregates and 
proposed oc 1/A°* 

For non-absorbing particles, van de Hulst [I] recom- 
mends the ratio of scattered intensity to equivalent refractive 
index of the particle cloud. The latter is given by 



m 



mo + 



where V is the volume of one particle and mo is the refractive 
index of the surrounding medium. The ratio again yields 

7.2. Fractal agglomerates 

Small particles dispersed in a fluid are commonly known 
to aggregate into large randomly branched structures. This is 
typical behaviour of soot in flames, for example. The 
method above relied upon K <u a However, Ku and 



Shim [287] point out that since the scattering by particles 
in an agglomerate is proportional to the square of the 
number, this assumption may no longer be valid. 

Because of the random nature of the agglomerate it has 
recently become popular to treat them as fractals. In this, the 
number of individual spheres N inside a cluster with radius 
of gyration R t is written 



N = K f (R 9 fro)* 



where r 0 is the radius of the individual sphere and D t is the 
fractal dimension. It is dubious whether the agglomerate is a 
true fractal, as its properties are not independent of scale. 
This "self-similarity" is a fundamental property of fractals 
which does not appear to be satisfied here. Nonetheless, the 
fractal description is proving to be a very convenient and 
powerful diagnostic tool. 

Quoted values of the fractal dimension vary. Depending 
upon the material it is usually in the range 1.5 to 2.0. Koylu 
et al. [288] point out that to fully characterize the agglom- 
erate, the constant K t must also be known. The specific form 
suggested for soot is [289] 



N = 2 



The authors claim that this is not sensitive to fuel type or 
flame condition. Sorensen and Feke [290] find Df = 1.8 but 
Kf = 1.7, and claim that the fractal morphology is constant 
for agglomerates ranging from assemblies of 10 up to 10* 
particles. 

The fractal dimension appears to depend upon the history 
of formation of the aggregate and can yield information on 
the formation process. Colbeck and Wu [291] found 1.4 £ 
Dt s l .96 for carbonaceous smoke clusters depending upon 
the history, and D f = 1 .08 for magnesium smoke which had 
the form of chains. Also. Colbeck et al. [292] found that the 
fractal dimension of carbonaceous smoke increases due to 
processing in clouds in the presence of water vapour. It rises 
to Of « 2 to 2.5. and the clusters behave more like solid 
spheres for light scattering. 

The scattered intensity from an agglomerate may be 
obtained from 



liq) = hnM.^^jPiqa^ig) 



with (do7oXl)o and P{qa) being the Rayleigh forward scat* 
tering differential cross-section and the scattering form 
factor for the primary particles. M 2 is the second moment 
of the size distribution and k % aM^ D \ S{q) is the structure 
function and q = 2k sin(0/2). 

The fractal dimension can be measured directly by plot- 
ting the scattered signal as a function of angle. Provided q = 
2k sin(0/2) is sufficiently large, it is found that 

S[q) oc q-» 

Sorensen et al. [293] claim that the structure function S(q) 
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has this form provided qR t > I- However, Bonczyk and Hall 
[294] state that q a 3 is required and propose extending the 
range of angles that can be used by moving to the ultraviolet. 

When qR t < I, the so-called Guinier regime, Sorensen ct 
al. (295] proposed a means of obtaining N and d. Here 

S(q) «■ 1 " ^q 2 ^ 

and plotting S(q) against q 2 yields a straight-line plot of 
slope R\. Plotting the size from scattering extinction </?se> 
against ff f results in a set of straight lines of constant r 0 and 
N. A log-log plot of /?se against R t gives lines of constant 
slope D,/3. A network is then established from which the 
three parameters can be found. However, the authors state 
that the results are dependent on the assumed form for the 
correlation function describing the distribution of mass in 
the agglomerate. 

Khlebtsov and Melnikov [296J and Sorensen et al. 12971 
have explored the structure function in more detail. They 
find 



S{q) = 



I 



qR $ < 1 Rayleigh regime 

1 2 

1 - ~(qR $ ) q&% < 1 Guinier regime 



(«*.) 



qR % » I Power-law regime 



(Faeth and Koylu (298) give 
S(*)«exp(-[*K,]*/3j 

in the Guinier regime, which reduces to the form given 
above for m[qR$4ZU 

Measurements in the Guinier regime can yield the overall 
size of the agglomerate while the power-law regime 
contains information on the structure. However, complica- 
tions arise due to a number of factors, including finite par- 
ticle size, polydispersity multiple scattering, and the lack of 
self-similarity. 

To allow for finite size, Khlebtsov and Melnikov [296] 
proposed a modified structure factor of the form 



1 1)sin[(0 f - l)arct 



S(q) = 



where M is the cluster mass and 




0 f (D f +l) 



Sorenson et a). [297] take into account the variation of 
cluster size and suggest 

1 



5( 9 )= 1-3^) 



where 



/*n - A. ^0 + '/Pi) 



and 



n{N) being the distribution of the number of particles in the 
clusters. 

As implied above, the cluster ceases to behave like a 
fractal when scales are considered which are too small (of 
the size of only a few primary particles) or too large (of the 
order of the cluster size). At the lower (Rayleigh) end, the 
scattering by a few particles can be calculated directly. At 
the higher (power-law) end the problem is tackled by the 
introduction of a cut-off function to the cluster density- 
density correlation function. If this takes the form h{r/Q = 
exp[ - irttfl where £ is the characteristic cluster size 
related 10 J? g , then Sorenson et al. [297] find 





where ,F| is a confluent hypergeomctric function. Fig. 14(a) 
shows plots of the variation of scattered intensity with (qRJ 2 
and Rg. 14(b) illustrates the sensitivity to the cut-off func- 
tion assuming the form h - e" 1 ^. Fig. 14(c) shows some 
experimental results for soot in a methane/oxygen flame. 

The sensitivity to cut-off function has also been explored 
by Bold et al. [299] They concluded that the polar diagram 
is not sensitive to its form, although it does have significance 
for cross-sections. 

A description of light scattering by agglomerates of small 
particles in terms of their fractal structure has been given by 
Berry and Percival (300). They used a mean field theory 
with scalar analysis. Their results suggested that multiple 
scattering within the agglomerate was negligible for D f < 2. 
For D f > 2 it dominates. This was explored further by 
Nelson [301] who concluded that the mean field approach 
is valid for assemblies of up to 250 particles. Colbeck et al. 
[292, 3021 reported that the theory of Berry and Percival was 
in agreement with their experimental results for A < 2. 
However, for D f > 2 the mean field theory predicted an 
increase in Q«* which was not observed. The results were 
more in agreement with Mie theory for a sphere. Similar 
conclusions were reached by Khlebtsov 1205], who treated 
the fractal agglomerates as porous spheres in the anomalous 
diffraction approximation. 

The increase for scattering with D f < 2 has also been 
reported by Dobbins et al. [303]. They suggest that the 
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Fig. 14. Scatter diagrams for fractal aggregates: (a) influence of 
fractal dimension for a structure factor derived from a Gaussian cut- 
off correlation function; (b) influence of (3 in the cut-off function for 
a structure factor derived from a Gaussian cut-off correlation func- 
tion with Dj — 1.75; (c) influence of the width parameter of a zero- 
order log-normal size distribution function (ZOLD) for a structure 
factor derived from a Gaussian cut-off correlation function with 
Df - 1.73. Experimental results are shown for a pre-mixed 
methane/oxygen flame at height (A) above the burner of 12 mm 
<0) and 18 mm ( + >. (After Sorensen ct al. (297J.) 

specific extinction coefficient can be described by 

c pA y m 2 + 2 / 

4mc*K f 2 D < / 3D f Y^ W ~ 1 f 

An example of the increase in the extinction over that of a 
volume equivalent sphere is illustrated in Fig. 15. 
Because a fractal aggregate withD, < 2 may be a tenuous 
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Fig. 15. Specific extioctioa coefficient against number of primary 
particles, n, forXp = 0.024, m = 1 .55 - 10.78, D t - 1.75 and K f - 
9.0. The dashed line is the Mic theory prediction for the volume 
equivalent diameter of the aggregate. (After Dobbins et al. (303).) 



structure, it is reasonable to suppose that the incident wave 
may propagate with little disturbance. This concept led 
Farias et al. [160, 161] to develop a simplified model based 
on the RGD approximation. They claim that their model is 
accurate to within 10% provided that \m - \\ s I, the 
primary particles are within the Rayletgh limit on size and 
(m 2 - l)/(m 2 + 2) is not approximated. The method is also 
said to be satisfactory for some cases with D t > 2, even 
though the agglomerate is more compact. 

Singham and Bohren [304] have pointed to different prop- 
erties between illumination with unpolanzed and polarized 
light In the former case the scattering by the agglomerate is 
insensitive to interactive effects or fractal dimension. For 
polarized light there is sensitivity to both, and the element 
Mi4 of the scattering matrix is a strong indicator of the 
fractal dimension. The agglomeration introduces depolari- 
zation into the scattered light. According to Lu and Sorensen 
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Fig, 1 6. Depolarization ratio versus number of primary particles per 
aggregate for linear chains, m - XJb + i0.6 and x^** 0.2. (After Lu 
and Sorensen [303].) 
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[305], the depolarization ratio is given by 

-0.6 

but tends to saturate, as seen in Fig. 16. 

A further complication to the simple theory arises due to 
resonances. Shalaev ct al. [306] and Markel et al. [307] 
suggested that there exist resonances caused by coherent 
interaction between the particles in the agglomerate. 
These can cause considerable enhancement of the optical 
properties. Iskander et al. 1308] showed that there is also o 
resonance when the length of the agglomerate is of the order 
of the wavelength. At such a point the simple theory could 
be in error by as much as 50%. The same workers have 
developed empirical formulas to describe optical absorption 
by fractal aerosol agglomerates [3091- 

The agglomerates may be anisotropic due to non-spheri- 
cal shape. The measurement of shape factors has been 
explored by Wu and Colbeck [310]. They define the shape 
factor as 



where p is the particle density, po is the unit density, D vt is 
the volume equivalent diameter and D K is the aerodynamic 
diameter. It is demonstrated that 



and 



Neimark et al. [3111 suggest thai, to allow for shape, the 
fundamental definition of the fractal structure should be 
adapted. They propose 



N = a y (Lyw^) D * 



where Lis the length and W is the width. Evidently f or y = 
1/2 the length and width have equal influence and D y = O f , 
the usual fractal dimension. Similarly, for y = I length 
predominates so tot D Y = D L and for y = 0 then D y = 
D w . For a growing aggregate they define an affinity expo- 
nent H through 



- - (*)" 



W 0 

where W 0 and U are the original or reference values at any 
time in the growth. Then they show that 



1 + ~ 



i + /y 



2 " 2/7 

We note that for D L = D W ~ D u then H = 1 . For aggregates 
growing in a combustion system they find Of = 1 .75 and H — 
0.91. 

Any anisotropy, form or intrinsic, results in a rotation of 
the polarization of the scattered light relative to that inci- 
dent For randomly oriented particles there is partial depo- 



larization. This is measured in terms of a depolarization 
ratio. A scalar theory cannot calculate these effects and a 
full vector analysis is required. Chen et al. 1312) and Markel 
et al. [313] have tackled this problem. 

Wide-ranging reviews of the electromagnetic properties 
of fractal agglomerates have been given by Charalampopou- 
los [314], d'Alessioet al. [315], Faeth and Koylu [298] and 
Shalaev (316]. 

■ 

7.3. Polarization measurement 

Anisotropy in particles can be of two types: intrinsic or 
form. The first of these is a property of the material and is 
independent of shape. A sphere may exhibit intrinsic ani- 
sotropy, for example. Form anisotropy is a property of the 
shape and is associated with non-spherical particles. A 
feature of anisotropy. whether it be intrinsic or form, is 
that it will induce changes in the polarization state of the 
scattered light relative to that incident This property can be 
used to study the nature or extent of the anisotropy. 

It is common to refer to vertical (V subscript) and hori- 
zontal (H subscript) polarization. Usually the measurement 
plane is horizontal so that vertical polarization is perpendi- 
cular (<1> = 90°) and the other is parallel <<|> = 0°). Four 
possibilities exist: 

1. The incident polarization is vertical and so is that scat- 
tered, then 

2. The incident polarization is vertical but that scattered is 
horizontal, then 



sea 



'hv 



3. Hie incident polarization is horizontal and so is that 
scattered, then 

Acs ~ 'hh 

4. The incident polarization is horizontal but that scattered 
is vertical, then 

Aca = 'vh 



For isotropic spheres / HV = /vh = 0. It is generally true that 
/kv — 'vh — the reciprocity condition. 
The depolarization ratios are defined by 

Pv — Arv/Jvv. Ph = 'vh^hh 

Typically they are found to be of the order 1%. 

Measurements of depolarization ratio for flame soot (e.g. 
Haynes and Wagner 1317]) find pv of the order \% whereas 
a prolate spheroid assumption would suggest about 4% for 
long chains at a typical refractive index of m « t .56 + i0.56. 
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The interactive chain model, however, predicts pv ^ 1%. 
Thus while Haynes and Wagner concluded that the 
observed depolarization ratio was due to intrinsic anisotropy 
and did not display evidence of chain formation, such a 
conclusion would be dubious with the chain model. Also, 
Mttller-Dethlefs (318] did find evidence in a counterflow 
diffusion flame for an increase of pv with time. 

It must be emphasized, however, that depolarization 
cannot distinguish between form and intrinsic anisotropy. 
Some other shape-dependent variable must be employed 
to do this, such as streaming birefringence (e.g. Jones and 
Wong [319]). 

It can be shown [320] that in the Rayleigh region the ratio 
MKcn i& the same for spheres and spheroids of the same 
volume. Thus the method of size measurement using the 
ratio of scattered intensity to extinction is independent of 
shape. 

7.4. Disymmetry and polar diagrams 

As the particles become larger, forward scatter increases 
relative to other angles. A simple measure of this is the 
disymmetry, denned as the ratio of the scattered intensity 
at two symmetrical angles. Popularly 45° and 135° are 
chosen. This yields a monotonic function of size up to x ** 
2, or about 0.4 \im in the visible. 

More information can be obtained from the complete 
polar diagram measured for both polarizations. As the par- 
ticle size increases forward scatter dominates and the polar 
diagram develops a convoluted lobe structure. For particles 
larger than about 4 jun and visible light (x > 20) the first 
few lobes around 8 = 0 are well represented by Fraunhofer 
diffraction. According to this, the width of the forward lobe 
is given by 

sin0= 1.22A/Z) 

while its height is proportional to D*. Thus this peak 
becomes increasingly dominant and narrow. The number 
of lobes in the polar diagram is of the order x|m|. 

For a polydispersc cloud the lobes due to different sizes 
occur at different angles. This has the effect of washing out 
the structure and making it featureless, except close to the 
forward direction where the strong peak remains. 

The polar diagram can be used to measure particle size. 
However, for single large particles, the structure requires 
increasing angular resolution which ultimately limits the 
maximum size that can be determined. For panicle clouds 
structure is not a problem, but the polar diagram may be so 
featureless that it is difficult to obtain sufficient information 
to recover a unique size distribution. This is especially so in 
the presence of experimental error. Further, at most angles 
the scattering is very sensitive to refractive index, which is 
an extra problem if this parameter is unknown. 

There have been two traditional approaches to recovering 
the size distribution. In the first a collection of polar 
diagrams are generated from Mie theory by using a range 



of parameters in an assumed form of the size distribution 
function. Typical functions employed are the log-normal 
and Rosin -Rammler distributions. Brown and Wohletz 
[321] have proposed the use of the Weibull distribution 
which naturally arises in fragmenting particle systems. 
The Rosin -Rammler distribution is one form of the Weibull 
distribution. The alternative has been to employ a direct 
inversion scheme, such as that of Phillips and Twomey 
(see Section 9). 

New techniques are continually being explored. For 
example, Min and Gomez [322] have examined the applica- 
tion of the Fourier transform of the polar diagram to single 
particles. This results in a frequency and phase correspond- 
ing to the number and position of the lobes. In particular, the 
frequency is found to be a linear function of size for 
diameters above 10 ujn, with additional information being 
provided by the phase. The calculations suggested that accu- 
racy should be excellent in dilute sprays. If the refractive 
index is known, the size could be measured to within 1%. 
Alternatively, an uncertainty of 2% in refractive would 
result in a sizing error of about 3%. Counting rates up to 
30 kHz were possible. The authors confirmed their method 
on individual particles in the size range 74,75 to 74.95 tun 
using an angular range from 10 to 17°. This size range is 
very narrow, suggesting that further study is required. 

Glover et al. [323] demonstrate that if a slightly out-of- 
focus image is taken of a particle, a set of fringes is observed 
which are regular in the angular range 20 to 70°. These are 
part of the scattering pattern arising from interference 
between front-reflected and back-refracted rays. The fringe 
spacing is directly related to the size through 



where x is the size and m is the real refractive index. A 
pulsed laser was needed to freeze the fringes, which move 
with the particle. 

A very novel approach due to Ludlow and Everitt [324] 
reformulates the Mie theory and the RGD approximation in 
terms of Gegenbauer polynomials, which have the advan- 
tage of forming a complete orthogonal set In consequence 
the polar diagram can be readily transformed into a set of 
expansion coefficients via a simple integral over the enure 
angular range. A plot of these coefficients against order 
constitutes the Gegenbauer spectrum, which has a number 
of interesting properties. The spectrum is much simpler in 
appearance than the original polar diagram and has a distinct 
cut-off at the order given approximately by 

/ co = 2.1x + 0.l2 

where x is the size parameter. An example is shown in Fig. 
17, where the C\ are the coefficients of the Gegenbauer 
expansion, h is not clear how useful this method would be 
for a polydisperse system. 
Polar diagrams have the advantage of being relative. They 



BNSOOCtD: <XP 415162QA_I_> 



A.R. Jones/ Progress in Energy and Combustion Science 25 (1999) 1-53 



29 




0 30 60 90 120 150 180 
Angk(dcg-) 

(a) 

100 r 




1 




0 4 8 12 16 20 24 
Order 

00 

Fig. 1 7. (a) Angular scattering pattern and (b) Gegenbauer spectrum 
for x = 10 and m = 1 .2. (After Ludlow and Everiu [324].) 

will yield panicle size independently of concentration. 
However, if the concentration needs to be known, then an 
absolute measurement of scattered intensity is required. 
Because this is absolute the laser output should be moni- 
tored, and the system requires calibration. To do this 
particles with known scattering cross-section and concentra- 
tion are needed. A gas is ideal for this purpose. Its concen- 
tration may be calculated for given temperature and 
pressure, and the scattering cross-section is known for 
many gases. As an example, for oxygen = 1.05 x 
10" 31 m 2 . d'Alessio [320] quotes a number of examples. 



7.5. Photon correlation spectroscopy 

A drawback of the methods described above is that they 
rely for their interpretation on a knowledge of the refractive 
index. A method of sizing based upon mechanical properties 
would be independent of this. The velocity is such a prop- 
erty. 

Small particles (< 1 |un> suspended in a fluid exhibit 
random Brownian motion as a consequence of molecular 
bombardment The more massive, the particle the less 
significant this effect is. Thus measurement of the random 
motion can yield size. While a frequency shift is involved 
due to the Doppler effect, the treatment is essentially dealt 
with by elastic scattering. For this reason it is sometimes 
referred to as quasi-elastic light scattering <Q£LS). 

The consequence of the Doppler effect combined with the 
random motion of the small particles is a broadening of the 
spectrum of the probing light beam. The application of this 



diffusional broadening to soot in laminar diffusion flames 
has been described by Lhuissier et al. [325] 

The frequency shifts are small, and while the spectrum 
can be obtained with a high-resolution spectrometer, an 
alternative is to find the broadening indirectly via the result- 
ing degradation in coherence. To measure this the scattered 
signal is compared with itself at increasing delay times, r. 
For a single particle size the resulting correlation function 
decays exponentially, the rate of decay increasing as the 
linewidth increases — or the particles become smaller. This 
technique is called photon correlation spectroscopy (PCS). 
The specific form of the correlation function is 

Sir) = constant x exp{-2q 2 D { t) 

where D f is the diffusion coefficient 
kT 



and k is Boltzmann's constant. T is the temperature and /x is 
the viscosity of the fluid. Also 



4tt 



0 



In reality there will be a size distribution. The correlation 
function is then the integral over a series of exponentials. 
Inverting this to the size distribution is a difficult process, 
currently the subject of much discussion. 

PCS is an established technique for small particles 
suspended in a stationary liquid, and commercial instru- 
ments are available. When used in flowing fluids, however, 
two complications arise because of the motion. The first is 
turbulence, which adds an extra random motion to the parti- 
cles and further broadens the scattered spectrum. The 
second complication is due to the bulk motion of the gas 
which carries particles across the test space. The finite tran- 
sition time is inverted to a characteristic frequency. This 
effect is the Doppler ambiguity. 

A number of authors have explored PCS for measure- 
ments on soot in a laminar flow, which removes one of 
the variables. In this case the correlation function for mono- 
disperse panicles takes the form 



Sir) oc exp(-2v 2 D f T- vV/Wg) 



the second term in the bracket being due to the Doppler 
ambiguity. To minimize this it is necessary to keep v/Wb 
as small as possible, v being the velocity and W 0 the beam- 
width (i.e. WVv is the transit time). Scrivner et al. (3261 and 
Weill et al. [327) have used this method to obtain soot 
particle size in flames. The latter authors recommend that, 
for a beamwidth of W 0 = 250 u,m, 

D «* 10 nm: v < 4.5 m s _l 

D *■ 100 nm; v< 0.5 ms' 1 
Weber et al. [328] have discussed the influence of concen- 
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Fig. 1 8. Illustration of an experiment to measure extinction. Scattered light passes to O and transmitted light to F. A pinhole at F allows the light 
to pass to the detector. 



tration, velocity, particle size and detector aperture on corre- 
lation spectroscopy. In general, they conclude that for dilute 
aerosols with t < Wq/v PCS works well, in agreement with 
the above. However, the method fails for t s Wrfv and non- 
dilute systems. 

A modified approach to PCS related to photon counting 
has been proposed by Gonzalez et al. (329} Their technique 
is based on monitoring the statistical distribution of time 
intervals* between arriving photons. The authors claim that 
this method can be superior at very low light levels. 

If the particles are non-spherical an extra degree of free- 
dom is introduced as they may rotate. The correlation func- 
tion then has two components which decay at different rates. 
Rarity [330] was able to separate these by measuring the 
mutual coherence function between two detectors. In this 
way he could measure both axes of spheroids. Xic et al. 
[33 1 ] used electron microscopy to measure the size of the 
primary particles in an agglomerate and then used PCS to 
provide an equivalent diameter from which the chain length 
could then be established. 

A review of the use of PCS for particle sizing in combus- 
tion systems has been given by Bernard (332]. 

7.6. Spectral extinction 

Efficiency factors were reviewed in Fig. 4. For particle 
sizing the most significant feature is that the extinction effi- 
ciency rises monotonically for small x. A reasonable guide 
isthatjcjm - 1 1 < 2. This suggests measuring the extinction 
as a function of wavelength— the spectral extinction 
method. Practically, the upper limit is of the order of 
I u.m diameter. For larger particles essentially oscillates 
around a constant value. The mean extinction then varies as 
the square of the size. However, the oscillation would make 
the measurement non-unique, except for strongly absorbing 
particles where it is damped. 

Care must be exercised in measuring extinction. A sche- 
matic experiment is indicated in Fig. 18. Light scattered 
from the cloud arrives at the point O and transmitted light 
at F. If the true extinction is to be obtained, negligible 
amounts of scattered light must pass to the detector. It is 
the function of the stop to ensure that as so. We know that as 
particle size increases so the scattering becomes more 



strongly peaked in the forward direction. Thus the accepted 
angular range becomes ever smaller, as does the aperture 
03) in the stop. Also, if the transmission is very high, then 
the measurement may be of a small change in a large back- 
ground. 

If there are temperature gradients and fluctuations of the 
gas, the transmitted beam is continually in motion. The 
angular displacement of the beam is proportional to the 
temperature gradient and the path length through the flame 
(although positive and negative gradients may offset each 
other to some extent). The problem associated with this is 
that either the position of the aperture would have to con- 
tinually be changed, or a larger aperture be used which 
would allow in scattered light. 

One way of overcoming this difficulty is to use a lens to 
image the particle cloud onto the detector. In this way 
deflections are brought back onto the optic axis of the 
system, and always arrive at the same point. Unfortunately, 
this also concentrates any emitted light at the aperture. To 
reduce the significance of this, phase-sensitive detection 
may be used together with a narrow pass interference filter. 
Alternatively a high-power laser may be employed. 

Since Gen is very sensitive to refractive index, this para- 
meter needs to be known with reasonable accuracy. Further- 
more, the method is most useful for scatterers with 
refractive index not too different to one, otherwise the 
measurable size range becomes too restricted. According 
to Shifrin and Perelman [333] the recommended upper 
limit is m - l.S. If the particles are suspended in a material 
other than a gas their refractive index relative to the medium 
must be used. 

At high concentrations multiple scattering becomes 
important However, extinction relates to the loss from an 
infinitesimally small solid angle around the coUtmated inci- 
dent beam. The scattering into this is zero. In reality, where 
a finite aperture is used, there will be some scattering 
received which will be influenced by multiple scatter. This 
will be small provided the aperture is sufficiently smalt. 
Multiple scattering will not be significant to extinction 
measurements provided appropriate care is taken to select 
a suitable aperture. 

Where multiple scattering cannot be neglected, a correc- 
tion factor has been proposed by Szyrnanski 1334). Here, the 



BNSDCC1D: <XP 4151620A__L> 



AM. Jones / Progress in Energy and 

transmissivity has the form 

and values for X are provided from near forward radiative 
transfer theory. The method is found to agree well with 
experiment. 

Wallers [335] applied the spectral extinction method to 
condensing water droplets in steam turbines. A probe was 
developed with the light source and filters at one end and the 
detector at the other, and a slit between them through which 
the droplets passed. The whole probe assembly could be 
inserted into the turbine chamber. The data were inverted 
to recover the size distribution using the Phillips-Twomey 
method (see Section 9). 

A method allied to spectral extinction has been applied to 
obtain mean particle size and concentration in very dense 
media by Gougeon ct al. (336]. This is a double extinction 
technique and uses the ratio of the extinctions at wave- 
lengths of 0.6328 jim in the visible and 337 u,m in the far- 
infrared. Guidt et al. [337] have validated the technique by 
making measurements on particles of known size and 
concentration suspended in a highly viscous gel. 

Another feature of the wavelength variation of extinction 
has been exploited by Winklhofer and Plimon [338] to 
monitor hydrocarbon fuel in internal combustion engines. 
Scattering and absorption at two wavelengths is used to 
discriminate between gaseous fuel and liquid drops. An 
output wavelength of 3.39 u.m from a He-Ne laser was 
chosen because this coincides with an absorption band in 
the fuel vapour. In the visible the fuel vapour is almost 
transparent. 

A possible extension to inhomogeneous particles has been 
considered by Graff et al. [339]. By studying the modifica- 
tions to the extinction spectrum due to the presence of 
metallic cores in coated particles, they claim that it may 
be possible to measure a coating thickness to within I nm. 

7.7. Fraunhofer diffraction 

The Fraunhofer diffraction pattern is very sensitive to size 
in a simple way, and polydisperse samples can be analysed 
with good accuracy. Refractive index does not appear in the 
diffraction equation and comparisons with Mie theory 
confirm lack of sensitivity to this for sufficiently large par- 
ticles and small angles. Thus the problem of knowing refrac- 
tive index not does arise provided a small error in size is 
acceptable. For larger particles in a size distribution the 
refractive index begins to affect the higher angles measured. 
Zhang and Xu [340] predict that errors can be in excess of 
10% due to use of the incorrect refractive index. This is 
roughly in line with the predictions of Jones [173]. 

A further advantage of Fraunhofer diffraction is that the 
governing equation is extremely simple, making calculation 
straightforward. As a consequence of these advantages 
Fraunhofer diffraction has become a popular and powerful 
method. 
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The basic apparatus required is also simple and was illu- 
strated in Fig. Kb). A parallel beam of monochromatic light, 
conveniently from a laser, traverses the particle cloud. The 
transmitted and scattered light are collected by a lens and the 
intensity is observed in the focal plane. The diffraction 
pattern is independent of the position of the particle in the 
beam and only depends on the scattering angle. (The light 
distribution in this plane is the Fourier transform of that in 
the test space, where the particle cloud is located, and the 
lens is often referred to as the "Fourier transform lens".) 

The intensity distribution can be measured by traversing a 
photodetector across the focal plane. This is time-consum- 
ing and it is preferable to use a detector array. This may be 
either linear or, as used by Malvern Instruments, segments 
of concentric circles. However, as multiple detectors are 
used, calibration is needed to ensure that their relative 
response is known [341, 342]. 

The instrument response can be checked by using parti- 
cles of known size distribution. An alternative is the use of 
calibration reticles which have large numbers of discs of 
well-determined size photographically etched onto glass 
1342]. 

The use of diffraction in two- phase flows has been 
reviewed by Bachalo |343]. The properties of the surround- 
ing fluid can have significant effects, especially if there are 
density gradients which can cause deflections of the trans- 
mitted beam. 

Hirleman and Dellenback (344] proposed a method of 
overcoming mis, in which the scattered light passed through 
a filter composed of an array of birefringent elements that 
could be made transparent or opaque to a particular polar- 
ization by the Faraday effect. The array was switched 
sequentially to scan the diffraction pattern with rings of 
any required radii. Transmitted light was collected by a 
lens and passed to a single detector, so that calibration 
was not needed provided the response was linear. The centre 
of the diffraction pattern was determined independently by 
separating a fraction with a beam splitter and sending this to 
a position indicator. In this way, allowance could be made 
for deflections along the light path due, for example, to hot 
gases. 

A similar concept was used by Miles et al. [345], although 
here conditional sampling was employed. The diffraction 
pattern was only recorded when the position indicator 
showed that the beam was central. Tests performed in flow- 
ing hot gas indicated that very large errors could occur if 
conditional sampling was not used. For example, a 45 u*m 
particle was sized as 135 u^n without sampling and as 
43 fim with sampling. 

7.8. Mixtures of small and large particles 

The scattered intensity close to the forward direction 
increases strongly with particle size. At 0° Fraunhofer 
diffraction gives oc D 4 t though at larger angles it 
becomes more like D z . This means that the scattering is 
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Fig. 19. Illustration of the method of displaying both the angular 
and wavelength variation of scattered intensity on a CCD camera: 
<a> apparatus: (h) appearance on screen. (After Card and Jones 
1347 1) 



dominated by large particles; one 100 u.m particle is equiva- 
lent to 10 000 of size 1 tun. Further, small particles no 
longer satisfy the criteria for Fraunhofer diffraction. 
Although this limitation can be overcome with the use of 
Mte theory, there is still a lack of angular information. For 
these reasons diffraction-based instruments were limited at 
the lower end until quite recently. The measurement of 
smalt particles has been incorporated by combining back- 
scattering and diffraction (346]. 

A different approach was used by Card and Jones 
[347. 348). This made use of the facts that for large particles 
scattering is strongly dependent on angle but insensitive to 
wavelength and vice versa for small particles. A lamp with a 
continuous spectrum was used as the light source and an 
optical system devised which was capable of displaying 
both the angular and wavelength variation on a CCD 
camera. This is shown in Ftg. 19(a). Scattered light passes 
through the slit in A and the test space is imaged onto the 
CCD via a diffraction grating. Unscattered light passes 
through the pinhole to provide a measure of the extinction. 
The form of the image is as suggested in Fig. 19(b). 
Evidently, the system is capable of delivering the spectral 
extinction method, but the authors only used extinction to 
indicate concentration. 

Scatter diagrams for 03 urn and 90 urn particles are seen 
in Fig. 20(a). The expected dependencies on wavelength and 
angle are clearly demonstrated. Fig. 20(b) shows an example 
of the scattering by mixtures of these sizes. By comparing 
the measurements with theory, the volume fraction of the 
smaller particles could be determined to within less than 
1%. Use of the Phillips-Twomey inversion method (see 



Section 9) indicated that the bimodal size distribution 
could be recovered satisfactorily. 

7.9. Dense particle clouds 

We have seen previously the problems which arise when 
attempting to measure particle size in dense systems 
because of multiple scattering, and have explored some of 
the methods used to attempt to correct single-scattering type 
measurements for this. However, interest is developing in 
ways of using the multiply scattered light directly. An early 
proposal was. for example, that of Deepak and Green (241 1 
who predicted the scattering by a radiative transfer approach 
akin to the Monte Carlo method. In their technique direct 
transmission and forward scatter at some distance offset 
from the incident beam were measured. 

In L1DAR a laser pulse is transmitted to a target area in 
the atmosphere and the properties of that region are deter- 
mined from the reflected signal. The influence of multiple 
scatter on this have been examined by Benayahu et al. [349] 
who considered only double scattering, which they claim to 
be sufficient. Here the direct reflection and backscatter at 
some offset distance were observed. 

A promising concept has been discussed by Zaccanti et al. 
[350]. who examined the influence of multiple scattering on 
the shape of a transmitted light pulse. They conclude that a 
time-gated scanning-imaging system should be able to 
delect absorbing structures inside thick turbid media. 
However, picosecond pulses are required. The propagation 
of pulses in dense media has also been examined by Kuga et 
al. f3Sl]. They find that the speed of propagation of a pulse 
can be reduced by as much as 50% as a consequence of 
multiple scattering events. Other studies on pulse spreading 
and broadening include those of Prahl et al. (25 1], van de 
Hulst and Kattawar [252], Shifrin and Zolotov (253] and 
Brewster and Yamada [254]. 

An alternative approach is to penetrate the dense medium 
with a suitable probe. A recent version of this idea is the 
endoscopic grating velocimeter-granulometer proposed by 
Cartellicr |352]. 



8. Particle 



v in • 



&/. Measurement of scattered intensity 

This leads to the method of measuring the pulse height as 
individual particles traverse a light beam. Unfortunately, 
away from the forward direction it is very sensitive to 
refractive index. Furthermore, it is an absolute measurement 
which leads to other problems as we shall see later. 

The intensity curve may oscillate and. so. not produce a 
unique value of size. However, the scattered intensity is 
integrated over the aperture of the collecting lens. Also, if 
a noo- monochromatic light source is used, there is Integra* 
lion over wavelength. The combined effect of these is that 
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Fig. 20. Diagrams showing ihc variation of scattered intensity with both angle and wavelength. The shades are representative of the logarithm of 
intensity, with white being the highest and black the lowest, (a) Pure samples of particles with representative sizes of 0.3 u.m and 90 *un, (b) 
Mixtures of these two sizes. The experimental results are on the left and theoretical ones on the right for the mixtures indicated. (After Card and 
Jones (3471.) 



any oscillations tend to be smoothed out and a reasonably 
monotonic variation of intensity with size is achieved [353]. 

If the refractive index is known the instrument response 
can be calculated for spheres. Often, however, it is not 



known and scattering by non-spherical particles cannot 
presently be calculated with certainty for other than a few 
well-defined shapes. However, as mentioned earlier, 
forward scatter is least sensitive to refractive index and 
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Fig. 21. Comparisons of the scattering patterns of randomly oriented non-spherical particles with equal volume and equal area spheres: (a) 
cubes; (b) flat plates. (After Pollack and Cuzzi 1 175].) 



shape. For example. Fig. 21 is taken from ihe work of 
Pollack and Cuzzi [175] and compares spheres and non- 
spherical particles. The calculations for the cubes and flat 
plates were based on a semi-empirical theory which was 
found to provide reasonable agreement with experiments. 
The curves clearly demonstrate the advantage of operating 
close to the forward direction. Unfortunately, it is often 
necessary in practice to use a considerable aperture in 
order to collect sufficient scattered light. Much of the advan- 
tage is then lost due to the range of angles seen. 

While forward scatter confers some advantage, it also 
raises difficulties. The first is that the unscattered light 
must be blocked by a stop, which eliminates the use of a 
smalt range of angles around 0 = 0°. This is unfortunate 
because scattering is strongest close to this angle. Apart 
from the loss of intensity, the sensitivity to shape and refrac- 
tive index is increased. Use of a highly colli mated source, 
such as a laser, minimizes the required stop size but it 
remains finite. A second problem is that the forward direc- 
tion is more sensitive to stray light scattered off optical 
elements and flare from the light source which passes the 
stop. There is always a balance to be struck, therefore, 
between a stop small enough to eradicate these sources of 
noise but which cuts off too much scattered light and a larger 
one with the opposite consequences. 

A third problem with forward scatter is that the detector 
looks along the incident beam and the test space can be quite 
long. It is only limited by the depth of field of the optical 
system. Spatial resolution will be best at 90°, where the test 
volume is given by the intersection of the incident beam 
with the depth of field of the detector optics. 

Scattering by a single particle is proportional to the inci- 
dent light intensity and depends upon its position in a non- 
uniform field: such as the Gaussian form of a laser beam. For 
systems which rely on received intensity this is a problem 
since the size will apparently vary with position. Thus the 



path of a particle is an important unknown. This is the 
trajectory problem. 

Perhaps the most obvious way to tackle this difficulty is to 
make the incident beam uniform in intensity; the so-called 
"top-hat** profile (e.g. Grehan and Gouesbet [354]). For 
example, the laser may be expanded and only the central 
portion selected by an aperture. Alternatively, a neutral filter 
may be constructed with an inverse Gaussian density profile. 
A holographic method has been used to the same effect 
[354]. The problem with all such methods is that they 
throw away much of the available light A possible way 
around this problem for one dimension is suggested by 
Xie et al. 1355], The original Gaussian beam is divided 
into two halves. Both halves are then rotated to provide 
their minor images and they are recombined, as shown in 
Fig. 22. This yields a reasonable approximation to a one- 
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Fig. 22. Dlustration of how two half Gaussian curves may be added 
to provide a "top-hat" profile. The brokeai curves are Gaussiam and 
the solid line is their sum. 
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Fig. 23. Cross-section of a Gaussian laser beam. 

dimensional "top-hat" profile while retaining all of the 
light. 

Another approach is to mark the centre of the incident 
beam. Typically this is achieved in one of two ways: either 
using two beams with the same wavelength but orthogonal 
polarization or having different wavelengths [356, 357]. 
One of the beams is expanded and the other is aligned to 
coincide exactly with its centre. Two detectors are required 
with filters to select either the polarization state or colour, 
and signals are only accepted which arise from both beams. 

The methods outlined above have the disadvantage that 
the incident intensity is reduced so that the scattered light is 
weaker. This need not be a problem with large particles or 
high-power lasers, but could create problems in marginal 
situations. There can also be problems with panicles skirting 
the edge of the test space which are only partly illuminated. 
The latter can be avoided in the two-colour system by 
accepting only signals with a given ratio of transit times 
in the two beams, but this restricts the count. 

An alternative approach is to invert the light scattering 
data using known properties of the structure of the test 
space. A simple technique for this was described by Holve 
and Self [358]. The principle is illustrated in Fig. 23. It is 
based on the fact that the incident laser beam has a Gaussian 
cross-section, so that the illumination of a particle depends 
upon its trajectory. A small particle following path A will 
see a higher illumination than one following path B and, 
consequently, the scattered light will be more intense. 
Indeed, for the smallest particles the light scattered out 
from path B may be so weak as to be below the detection 
limit of the system and they will not be seen. A large particle 
following path B may scatter the same intensity as a small 
particle along path A. The biggest signals come from the 
largest particles following path A. 

In the inversion procedure the test space is divided into a 
number of zones within each of which the incident intensity 
is approximately constant. The maximum possible pulse 
height is given only by particles in the largest size class 
passing through the central zone of maximum brightness. 



The next class of pulse heights is provided by the largest 
particles passing through the adjacent zone, which is one 
step down in brightness, and the next particle size class in 
the central zone, and so on. In this way the number of pulses 
of any height is related to the intensities in the test space 
zones by a triangular scattering matrix. It is assumed that the 
flux of particles is uniform, such that the size distribution 
and probability of arrival are independent of the path. Then 
the matrix can readily be inverted to retrieve the numbers of 
particles in each size class provided the incident intensities 
are known; that is, the structure of the test space. This can be 
obtained by calibration with monodisperse panicles. More 
recently the matrix inversion has been replaced by an inte- 
gral deconvolution (359], and the structure of the test space 
calculated from a knowledge of the laser beam and apertures 
in the apparatus. 

In theory there is no upper limit on the particle size. In 
practice there may be limitations due to finite beam width 
and apertures, and non-linear detector response or satura- 
tion. Ultimately it becomes more practical to resort to 
simpler and more straightforward methods, such as photo- 
graphy, holography or schlieren techniques (e.g. Birch (360] 
and Harada and Murakami (361]). 

Ultimately the smallest measurable size will be limited by 
noise due to scattering by gas molecules. It is worth 
mentioning in this context that Ishikawa et al. (362) claim 
to be able to detect nano-size particles in water by employ- 
ing photon counting. Their method compares moving 
averages of the statistically fluctuating signal. A particle is 
present when the average rises above a threshold set by the 
average of the background. They were able to size particles 
in the range 0.067 u.m to 0.2 \im. 

The operation of particle counters requires the presence 
of only one particle in the test space at a time. This limits the 
upper concentration to roughly the reciprocal of the volume 
of the test space. For example, a spherical volume with 
diameter I mm could accommodate about 10 9 particles 
per cubic metre. There is. in principle, no lower limit to 
the concentration but the counting time may become too 
extended. The particles arrive at random and are subject to 
Poisson statistics according to which the standard deviation 
is proportional to the square root of the number of counts. To 
determine the number of panicles of any one size class to 
within 1% would require that 10 000 were detected. 

The size distribution is built up particle by particle so 
there is no inversion problem. However, a temporal average 
results. If the spatial average is required this may be recov- 
ered if the velocities of the particles are known. Holve and 
Davis (359] used the known structure of the test space, and 
measured the time taken for the signal height to change 
between two fixed fractions of the maximum. Since this is 
a known distance, the velocity is found. Alternatively, laser 
Doppler may be employed. 

In a combustion system problems arise due to beam fluc- 
tuations caused by refractive index gradients in the gas. To 
avoid this a gating method may be employed in which both 
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Fig. 24. Schematic diagram of apparatus to discriminate between fibrous and non-fibrous particles. (Alter At-Chalabi ct al. [367),) 



the laser beam centre line and the centre of the sample 
volume are constrained to coincide before a signal is 
accepted. 

An attempt was made to minimize the influence of non- 
spherical particles by Hardalupas et al. [363] These authors 
examined the diffraction patterns of individual particles and 
rejected any signals where the pattern had a large aspect 

ratio. 

Finally, the measurement of intensity is an absolute 
method. It is sensitive to changes in incident intensity, due 
to variation in source output, unclean optical elements or 
beam blocking by particles along an extended path through a 
cloud. The scattered light may suffer attenuation for the 
same reasons. 

There is also the need for calibration using particles of 
known size. Ideally the particles should be spherical of well- 
defined size. Unfortunately, these are not always readily 
available at sizes of a few micrometres. Polystyrene latex 
spheres are popular, but these come in aqueous suspension. 

Sample ink 



To be used in air the suspension must be diluted, sprayed 
and evaporated to leave behind individual particles. 

8.2. Non-spherical particles 

Recently, attempts have been made to extend the use of 
particle counting to the characterization on non-spherical 
particles. Their influence may be detected from the structure 
of the scattered pulse, the lack of azimuthal symmetry or the 
polarization states of the scattered light. 

The variation in response due to a large particle crossing a 
smaller beam has been used to provide information on the 
structure of the particle surface (e.g. Buckle and Bushnell- 
Watson [364] and Ohba and Isoda 1365]). As the beam 
strikes different facets on the surface the scattered light 
fluctuates in intensity. A similar concept is the so-called 
"slit scan*' method [366] in which a large particle is 
scanned by a narrow beam to reveal structure. 

Of particular relevance to industrial safety, especially 
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with regard to asbestos and mineral wool, is the detection of 
fibrous materials in the atmosphere. Examples include the 
polarization method of Al-Chalabi et al. (367] and the aniso- 
tropic azimuthal scatter method of Kaye et al. [3681 In the 
former, Light from a laser is converted to circular polariza- 
tion. Light scattered from a spherical particle is similarly 
circularly polarized, but from a fibre it is elliptically polar- 
ized. This change is detected by passing the scattered light 
through a quarter-wave plate and then to two photomulti- 
pliers which observe orthogonal polarizations. An increase 
in signal at polarization orthogonal to that incident signifies 
the presence of a fibre. The method is illustrated in Fig. 24. 
Rood et al. (369J used a similar concept by observing the 
orthogonal polarized intensities, but from panicles sampled 
through an aligning nozzle onto a slide. 

The instrument designed by Kaye et al. 1368] is shown in 
Fig. 25. Here particles are passed through one focus of an 
ellipsoidal cavity. They are illuminated by a laser beam 
which enters through a small aperture and is reflected onto 
the particles by the mirror M. Scattering is observed over a 
wide range of angles at the second focus. The presence of 
form anisoiropy is indicated by azimuthal asymmetry in the 
scattering partem. Originally, three detectors at the vertices 
of an equilateral triangle were used. A CCD camera has also 
been employed with the promise of the ability to character- 
ize quite complex shapes [370, 3711. More recently, 
attempts have been made to compare the observed patterns 
with theoretical predictions using the RGD approximation 
[372, 373]. The agreement is encouraging. 

Using the original three detectors, Kaye et al. [374] have 
developed an automated method capable of monitoring up 
to 10 000 particles per second. The shape factor is obtained 
from 



40.81 



where the three values of Ej are the measurements at each of 
the detectors and £ is the average. The mysterious value of 
40.81 is to ensure that the maximum value of A f — 100. Of 
course for spheres, A f = 0. 
The same instrument can be used to obtain size from 

3 I 

where £* is the forward-scattered signal This measurement 
can be used to obtain the spherical equivalent size for refrac- 
tive indices in the range 1.33 ^ m ^ 1.38. 

Similar ideas have been explored by other authors. For 
example, Sachweh et al. [375] also examined the asymmetry 
of azimuthal scattering, though at specific angles. Marshall 
[376] studied the diffraction pattern passing through 
specially designed filters to the same effect. In a slightly 
different approach, a fibre detector was developed by Brenci 
et al. |377] based upon measuring extinction and scattering 
at two angles. 



Other examples of the use of polarization include the 
study by Belanger et al. [378], who illuminated a mixture 
of glass fibres and spheres designed to simulate skin or an 
injection-moulded polymer blend. By using polarized inci- 
dent light and taking an image of the particles through a 
crossed polarizer, distinct patterns were observed which 
were indicative of the structure. Sierra and Mora [379] 
and Kim and Martin [380] studied the depolarization due 
to spheroids, in the latter case as a function of wavelength to 
develop a diagnostic method for interstellar grains. 
Mishchenko and Hovenier [381] performed a more general 
analysis for randomly aligned non-spherical particles. 

A further application to non- spherical particles was 
described by Card and Jones [382J. This method used 
cross-polarization techniques to discriminate between coal 
and fly ash. The scattered light from the particles was 
collected at an angle close to the backward direction and 
analysed for polarization state using two photomultipliers in 
the same manner as in Fig. 24. By plotting the cross- polar- 
ized intensity against particle size, it was found that the two 
types of particle were readily discriminated. More recently 
[383] it has been demonstrated that there is a good correla- 
tion between the ratio of scattered intensities for the two 
orthogonal polarizations and the carbon content of the 
combusted coal particles, as suggested by Fig. 26. 

Umhauer and Bottlinger [384] and Sasse (385] have made 
observations on single irregular particles suspended in an 
electromagnetic balance. The former authors compared 
measurements with calculations obtained by applying 
geometrical optics to slices of the particle in turn. They 
found thai the ratio of the width to height of the scattered 
spectra was constant. This leads to possible means of elim- 
inating the effects of shape and structure. Sasse [385] and 
Sasse et al. [386] have examined the phase function and 
cross-polarization from absorbing particles and means of 
directly measuring the albedo. 

Another method of suspending individual particles was 
propounded by Doornbos et al. [387] This uses the radiation 
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Fig. 27. Illustration of the coordinates for the Doppier shift equa- 
tion. 

pressure from two laser beams to trap single biological cells 
fox optical study. The beams came from two 100 mW diode 
lasers with output wavelength 780 nm. Particles up to the 
order of 7 u,m could be supported. 

An apparatus for measuring all the components of the 
scattering matrix has been described by Kuik et al. |388] 
A number of studies have investigated the most suitable 
elements of this matrix for characterizing non-spherical 
particles. For example. MLshchenko et al. (91) applied the 
r-matrix method to the problem of two spheres in contact 
They concluded that the only distinct manifestations of non- 
sphericity were that (M 22 fMn) * I, (Afjj/Af,,) * {MJM U ) 
and non-zero backscattering depolarization ratio. Martin 
(389) studied slightly non-spherical panicles using a pertur- 
bation method and obtained solutions for the matrix 
elements as expansions in spherical harmonics. In particular 
he found that for a rotationally symmetric particle 

"» " = tan 21* 

where \<t> is the angle between the azimuthal angle of the 
particle axis and the direction of the incident polarization. 




<b) 

Fig. 28. Determination of the beat frequency in IDA from inter- 
ference. Id (•) the scattered frequency is a function of angle; in (b) it 
is independent of angle. 



8.3. Laser Doppier methods 

The need to measure velocity in particle counting 
suggests the use of optical anemometry. A suitable vehicle 
for this is the Doppier shift suffered by light deflected by a 
moving object. In the coordinates of Fig. 27. the scattered 
frequency seen by a stationary observer at the angle 0 is 

f*i = »W| 1 + — sin ~ stnf a - - J 

where c is the speed of light. This equation has two implica- 
tions. First, since u «; c, the change in frequency is very 
small compared with that of light. Secondly, the scattered 
frequency varies with angle. 
The maximum possible change in frequency is 

K» * c 

Thus for a velocity of 10ms" 1 , say, we find 

— = 7x I0~ 8 
v 

Such a small change would be difficult to detect directly. 
However, for a typical optical frequency of 5 x 10" Hz, 
then Av = 3.5 X 10* Hz. This is coming into a range 
which can be readily be detected electronically, and it can 
be reduced further by appropriate choice of scattering angle. 
It can be observed by mixing the scattered wave with 
unscattered light to obtain the beat frequency from their 
interference, as in Fig. 28<a). The beat frequency varies 
with angle. Consequently, uncertainty in angle leads to 
error, and the received frequency can vary across an aper- 
ture. Also, it is often required to make measurements at or 
near the forward direction (0 = 0°) where the beat frequency 
becomes zero. 

These problems can be avoided by the arrangement in 
Fig. 28(b). Here the particle crosses the interference pattern 
produced by the intersection of two light beams. Theoreti- 
cally this can be treated in two ways. Either the Doppier 
shift can be calculated for the wave scattered from each 
beam independently and the two waves added to give the 
beat frequency, or it may simply be considered in terms of a 
particle crossing a fringe pattern. In the latter case it is easily 
seen that the received frequency is given by 

v = w/A r 

* 

where A f is the fringe spacing. The beat frequency treatment 
yields exactly the same answer. In this crossed beam 
arrangement the received frequency is independent of scat- 
tering angle. Also, the fringe spacing can be measured 
directly, reducing the error due to angular measurement. 
The test space is small, being given by the volume of inter- 
section of the two beams, so that spatial resolution is 
improved. 

A consequence of the size of the test space is the maxi- 
mum allowable concentration. There mast be only one par- 
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ucle at a time in this volume otherwise interference occurs 
between the scattered signals. For focused laser beams a 
diameter of 100 pm is reasonable. Then the maximum 
concentration is of the order 10 12 nT 3 . This may be reduced 
in forward scatter depending upon the achievable depth of 
held. 

The general form of the scattered signal is 

= A + B cos(2m* + <f>) 

The frequency v yields the velocity of the particle. For this 
reason the technique is referred to as laser Doppler anem- 
ometry (LDA). This subject has been reviewed by Durst et 
al. {390] The amplitudes A and B and the phase are functions 
of the particle size. Particle sizing techniques have been 
developed based on the measurement of A, the visibility 
{BIA) and cV 

The most obvious way to measure the size is via the mean 
intensity A. Such an instrument is a particle counter of the 
type described in Section 8. 1 with added velocity measure- 
ment and the advantage of good spatial resolution. Apart 
from this it suffers from the same problems as the other 
absolute techniques. The method has been explored, for 
example, by Yule el al. [39 1], Ereaut [392], Orfanoudakis 
and Taylor (393] and Tcrhaag and Renz [394] using either 
pointer beams or top-hat beam profiles [354]. The combina- 
tion with two-colour pyrometry to also obtain temperature 
has been discussed by Finke et al. [395]. 

The possible advantage of using the visibility or the phase 
is that they are relative measurements. They are independent 
of the instrument provided its response is linear. 

8.3.1. Signal visibility 

The visibility as a function of particle size can be calcu- 
lated from Mie theory (e.g. Jones [396] and Eliasson and 
Dandliker [397]) or using common approximations (e.g. 
Chu and Robinson [398], Bachalo [399] and Pendleton 
[400]). A very simple result is that for large particles and 
collection apertures the visibility in the forward direction 
id = 0°) reduces to 

27,(770/^) 

where Mz) is a Bessel function and A f is the fringe spacing 
in the interference pattern [401-403]. 

There are several major weaknesses of this method. The 
equation above shows that the visibility is initially a mono- 
tonic function decreasing with size, but then it oscillates. 
Thus the dynamic range is limited since beyond some 
point (typically D *• Xf) the size is not uniquely determined. 
Detailed Mie calculations (e.g. Hong and Jones [404]) also 
reveal extreme sensitivity to refractive index at even quite 
small angles away from the forward direction. If this prop- 
erty is not known, then very large errors can ensue. In 
addition, experiments have shown that the visibility is too 
shape -sensitive for the method to be employed on anything 
other than near spheres [405]. 



39 

For non-absorbing particles the variation of visibility with 
particle properties and angle increases with the scattering 
angle until, near the backward direction, measurement is 
almost impossible [406]. On the other hand, for absorbing 
particles the visibility becomes unity for all angles greater 
than about 40° and there is no information. Since in many 
practical situations access is only available near the back- 
ward direction, these are important limitations. 

Negus and Drain [407] suggested ways to overcome some 
of these problems from a study of the effects of aperture type 
and size and polarization state. In backscatter they were able 
to obtain a monotonic relationship in the range 0 to 250 urn. 
Measurement of backscatter was further explored by Azar 
and Verurice [408], who produced maps of the real part 
versus the imaginary part of refractive index to delineate 
regions in which the visibility method is usable for incident 
circular polarization. They investigated a range of fringe 
spacings. aperture sizes and off-axis angles. 

The visibility can never be greater than the fringe contrast 
in the interference pattern. This can be degraded due to 
unbalanced incident beam intensities. This can be caused 
by blocking of one beam or unequal blocking on passing 
through an extended particulate medium. Lack of temporal 
coherence can also result in the same effect These phenom- 
ena and their significance have been explored, for example, 
by Yianneskis |409], Kliafas et al. |410, 41 1] and Lagranja 
et a). [412]. 

Also of significance may be refractive index gradients. 
These introduce two effects. First, phase shifts in the inci- 
dent beams cause the fringes to move within the test space. 
Secondly, deflections of the beams cause the test space to 
wander. In anemometry these lead to false signals and errors 
in the measurement of velocity. An apparent velocity can be 
recorded from a stationary particle [413]. A similar problem 
has been reported by Ancimer and Fraser [4 J 4], who also 
suggest means for minimizing the error. 

An interesting alternative to fringes of varying amplitude 
is fringes of varying polarization. Avery and Jones [415] 
showed that if the laser Doppler arrangement was modified 
by passing one beam through a quarter-wave plate and the 
other through a three-quarter plate, the result was a contin- 
ual rotation of the plane of polarization across the test space. 
The period of the rotation was the same as the fringe 
spacing. Massdi et al. [416, 4171 have used this method to 
measure the ratio of intensities scattered for two orthogonal 
incident polarizations. They claim to be able to measure 
drop sizes in sprays of up to several hundred micrometres. 

8.3.2. Signal phase 

A technique which overcomes many of the problems 
discussed above is to measure the phase of the scattered 
signal. This was proposed by Durst and Zare [418] and 
came into prominence following the work of Bachalo and 
Houser [419] and Saffman et al. [420]. A history of the 
development of phase Doppler anemometer has been 
given by Hirieman [421]. 
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Test space 




Fig. 29. Schematic diagram of phase Doppler particle analyser. 
(After Bachalo |424].) 

As the particle moves in the lest space formed by the two 
laser beams, the scattered light generates a set of inter- 
ference fringes which move across the face of any detector. 
The resulting signal at a fixed point in space oscillates in 
time as the fringes move past. Because of the way in which 
the fringe separation varies in space, the phase <£ varies 
similarly. This change of phase can be observed by compar- 
ing the signals received by two or more spatially separated 
detectors. 

The rate of change of phase with angle generally 
increases with particle size. Under appropriate conditions 
it is found that the change is linear with sufficient slope to 
provide suitable accuracy for particle sizing. It is also linear 
over a wide range and the dynamic range can be large — 
typically 0.5-3000 u.m. The accuracy at the smaller sizes 
has been discussed by Sankar et al. [422]. who found that for 
O < 10 nm an accuracy of ± 0.3 ujn is possible. Further 
extension to sub-micrometre particles is considered by 
Naqwi et al. |423]. 

There is a limitation, however, because once the phase 
becomes 2tt the wave looks identical to one with a phase of 
zero. This can partially be overcome by the use of multiple 
detectors. In the instruments currently available up to three 
are used, typically as shown in Fig. 29 [424]. The scattering 
angle is either 30° or 150° out of the plane of the laser 
beams. The phase difference between the two outer detec- 
tors is greater than between either of these with the central 
one. It can be arranged that while the outer phase difference 
varies from zero to 67r„ the centra) change is up to 27r. The 
outer detectors provide good sensitivity to size, while the 
central detector indicates which of the three cycles is being 
observed. 

A generalized treatment of light scattering out of inter- 
ferometric devices has been published by Naqwi and Durst 
[425]. This allows for arbitrary polarization, non-plane 
waves and Gaussian beams, and for integration over a vari- 
ety of aperture shapes. Certain polarizations may assist in 
the production of a linear response by suppressing reflected 
light. Further analysis by Naqwi et al. [426) has suggested 
angles and apertures suitable for the measurements of metal* 
lie particles as small as 20 nm. transparent particles in the 



sub-micrometre range and particles with refractive index 
close to that of the surrounding medium. 

The phase Doppler method has the advantage of being 
intensity-independent provided there is sufficient visibility 
to yield a measurable signal. It will not suffer from beam 
blocking nor need particles be restricted to the central zone 
of high fringe contrast. However, when particle counting, a 
correction is necessary to allow for the bias towards larger 
particles because of their larger scattered signal (409). 

In the phase Doppler method the fringe spacing can be 
less than the particle size and the test volume can be kept 
small. Use of the off-axis angle of 30° also makes the inter- 
section area smaller. Bachalo [427] has reported that 
concentrations of up to 10 12 m" 3 can be used. 

The trajectory of the particles through the test space can 
lead to errors because the sample volume is a function of the 
particle size, as in all particle counting techniques with non- 
uniform intensity profiles. Schaub et al. [428) employed an 
arbitrary beam theory to explore this and found that phases 
are smaller for trajectories through the edge of the test space 
away from the receiver, and vice versa. This leads to an 
artificial apparent broadening of the size distribution. Scat- 
tering resonances were found to have a similar effect. The 
phase errors become particularly significant when the beam- 
width is not large compared with the particle diameter, as 
discussed by Albrecht et al. [429, 430). The signals are only 
valid for the short period when the beam lies along the 
panicle diameter. 

Xu and Tropea [4311 have proposed alignment proce- 
dures to minimize such errors arising from the Gaussian 
beam profile. Another suggestion for avoiding trajectory 
problems when the diameter of the test space is similar to 
that of the particle is given by Hardalupas and Taylor (432). 
This involves the application of validation criteria based 
upon rejecting signals for which the phase lies outside the 
expected range in the three-detector system. 

An alternative which moves away from the phase Doppler 
concept has been proposed by Borys et al. (433). In this a 
frequency shift is added to one of the laser beams and the 
Doppler burst moves as the particle travels across the test 
space. This movement is simply related to the size. 

Means of determining particle concentration and flux 
have been discussed, for example, by Albrecht et al. 
[429. 430). Qin and Sorrunerfeld [434] and SommerfeW 
and Qin [4351- The difficulty here is that, while the numbers 
of particles can be counted, there is some uncertainty in the 
cross-sectional area through which they pass because it is 
size-dependent. This arises because, as mentioned 
previously, a large particle near the beam edge will scatter 
a measurable signal whereas a small one will not. The 
interpretation is that the larger particles are falling through 
a larger area. Thus the size and trajectory are also important 
here. 

Tropea et al. [436) have attempted to resolve this diffi- 
culty with the design of a dual-mode phase Doppler anem- 
ometer. The principle is illustrated in Fig. 30. Two pairs of 
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Fig. 30. The optica] geometry of the dual -mode PDA. (After Tropea 
et al. [436U 

laser beams converge onto the lest area in orthogonal planes 
so that two sets of interference fringes are generated at right 
angles to each other. The two sets could be discriminated, in 
principle, by using different wavelengths or polarizations. 
The detection system resides at 30* out of the scattering 
plane of one set of beams, as in the standard instrument. 
However, an additional pair of detectors measures the 
component of a particle's velocity orthogonal to the stan- 
dard system, thus providing complete information on its 
trajectory. The effective width of the measurement volume 
is ascertained from the length of the Doppler burst. 

Off-axis measurements raise the questions of shape and 
refractive index. It would be expected that the method is 
sensitive to these. Currently the particles must be spherical 
or nearly so. and the refractive index may need to be known, 
in order that the instrument response can be calibrated or 
calculated. Schneider and Hirieman [437] have also pointed 
out that there may be significant errors in phase due to the 
presence of refractive index gradients, although the response 
remains linear. Nonetheless, the inhomogencity results in a 
blurring of the recovered size distribution and Koser and 
Wriedt [438] have proposed an iterative technique for 
removing this effect. 

The influence of shape has been investigated by Doicu et 
al. (4391 using the extended boundary method to solve the 
problem for spheroids. They find that an eccentricity of 0.05 
causes a phase error of 5%. 

In two-phase flows small tracer panicles are often used to 
indicate the flow of the gas. An example would be the flow 
of the gas within and around a liquid spray. It would be 
convenient to measure both particle sizes with the same 
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optical system but, unfortunately, the scattered signal from 
the large drops usually swamps that from the seed particles. 
Friedman and Renksizbulut [440] have avoided this 
problem by adding an absorbing dye to the large drops to 
reduce their scattering cross-section. 

8.4. Size from the measurement of inertia 

In a diverging flow a particle accelerates at a rate depen- 
dent upon its inertia. Injecting a particle into a tailored flow 
pattern and measuring its velocity will yield its size. There 
are a number of instruments based upon this principle. Such 
a measurement is purely mechanical and does not depend 
upon a knowledge of the refractive index. The size recorded 
is that which determines the mechanics of the interaction 
between the particle and the carrying fluid. This is the 
hydraulic or aerodynamic diameter. In many applications 
in fluid flow this is the required parameter. 

An approach proposed by Tate et al. [441 J uses a fibre 
optic probe inserted into the flow. The probe has a hemi- 
spherical end which acts as a bluff body causing divergence 
of the fluid and having a stagnation point a short distance 
ahead of it where the fluid velocity is zero. Due to its inertia, a 
particle approaching the probe will not have zero velocity at 
this point. The more massive the particle, the closer its velo- 
city will be to that at a distance remote from the probe. Hence 
measurement of the velocity at this point yields the mass and, 
consequently, the size provided the density is known. 

Light is transmitted along the fibre and focused onto the 
stagnation zone. Reflection from the particle is collected and 
passed back to a detector where it is mixed with unscattered 
light. A beat frequency is then observed arising from the 
Doppler shift. 

The minimum measurable particle size is governed by 
two factors. There must be sufficient light intensity to 
provide a good signal. Given this, the velocity at the stagna- 
tion point musl be sufficiently different from zero, which 
suggests a minimum inertia. Similarly there is a maximum 
inertia because the velocity must be different from that 
remote from the probe. Tate et al. (44 1 ] claim that a measur- 
able size range in the ratio of ten to one is feasible. 

There is some error due to the detection of particles not at 
the stagnation point arising from the finite depth of field. The 
upper concentration thai can be used is limited by multiple 
particles in the effective test space. 

In turbulent flow there are orthogonal velocity compo- 
nents. This means that a measured particle will not generally 
be travelling directly towards the probe. Only a component 
of its velocity will be obtained and it will be judged to be too 
small. However, results due to Vincent et al. [442] imply 
that this effect will not be serious in the range of Reynolds 
number between 3300 and 10*. 

An interesting technique based upon aerodynamic proper- 
ties has been suggested by Roll et al. [443]. In this particles 
are optically levitated on a laser beam and their behaviour 
monitored when the beam is momentarily interrupted. The 
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Fig. 3). Use of total internal reflection for the study of particle 
deposition. The refractive index of the half cylinder (n 2 ) is greater 
than that in the chamber (n,) and the angle of incidence (6) is greater 
than the critical angle for total internal reflection. (After Polverari 
and van de Ven |464].) 

method is not restricted to spheres. The authors used it to 
study the evaporation of single and two-component drops in 
the size range 10-40 jim. 

8.5. Particles on surfaces 

A new area which has sparked much interest is scattering 
by particles on or near surfaces. This has arisen primarily 
because of the problem of contamination of surfaces in 
micro-electronic circuit production. There is a need to detect 
and identify panicle contaminants. Recent papers in this 
area are by Lindell et a). [444]. Muinonen et al. [445], 
Videen [446. 447). Johnson [448, 449), Jakeman [450] 
and Saiz et al. [451]. 

Videen et al. (452] have compared four different theories. 
They all agree in the far field, but approximations break 
down when the spherc-to-surface distance is less than the 
wavelength. More recently Videen and co-workers have 
published a series of papers where they develop a model 
which only requires the expansion coefficients for the scat- 
tered wave of an isolated panicle. The particle can be of any 
shape. The theory treats the particle and its image subjected 
to three fields, the incident and reflected fields and a term 
representing the interaction due to the reflection of scattered 
light (453-455]. 

Typical further techniques include direct backscatter plus 
scattering of reflection from a surface [456], the discrete 
sources model |457J, ray tracing [458]. image theory 
[459] and multipole expansion f I9J. Some special applica- 
tions which have been discussed include anisotropic par- 
ticles [460]. cylinders (455]. the effects of particle (tensity 
|461] and sizing metallic particles [462]. While the above 
have been concerned with metallic surfaces, Lis with et al. 
(463) have specifically examined scattering by small 
spheres on silicon and oxidized silicon substrates. 

In an interesting variation. Polverari and van de Ven (464] 
studied the deposition of panicles on a transparent surface 
using frustrated total internal reflection. The incident wave 
was inside the material, the angle being sufficient for tout 



internal reflection. Panicles were deposited on the outside 
and scattered the evanescent wave which penetrated the 
surface. The principle is seen in Fig. 31 . Since the penetra- 
tion of the evanescent wave depends upon the angle of 
incidence, variation of this angle allowed different levels 
in the panicle layer to be observed. 

9. Direct inversion to size distribution 

To obtain the size distribution for scattering from a cloud 
of particles, it is necessary to inven the light scattering data. 
There are two approaches. The first is the "best-fit'* method 
in which the scattering properties are predicted by using 
presumed size distributions and compared with measure- 
ments to find the closest comparison. Commonly, simple 
two-parameter distributions are assumed such as Rosin- 
Rammler and tog-normal functions. The problem then 
reduces to the determination of only two numbers. This 
works well provided that the true distribution is similar to 
that presumed. However it often breaks down, especially if 
the true distribution is multi-modal. 

The alternative is direct inversion, in which an integral of 
the form 

fix) = JgU.y)n(y) dy + e{x) 

is solved to obtain n(y).f(x) is the measured function and 
e(x) is due to errors of measurement gix 9 v) is the kernel of 
the integral and is a known function which may be calcu- 
lated from an appropriate theory. 
Specific forms of the terms in the integral equation are: 

1. Fraunhofer diffraction 
fix) = 

. 4 r2/,(**0)"p 

<v) = n(D} 
where D = 2R. 

2. Spectral extinction 
fix) = tf MI (A) 

irDr 
n(y) = n(D) 

3. Photon correlation spectroscopy 
fix) = Sir) 

gU.y) = exp(-2* 2 0 f T) 
My) = i»(D f ) 
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Rg, 32. Ray diagram illustrating the principle of the generation of 
two pulses by a spherical particle passing through a narrow laser 
beam. (After Onofri ct al. {505, 5061.) 



En this case the inversion takes the form of a Laplace 
transform. 

We may write the integral equation in matrix form as 

In principle, if there is little or no error (£, — 0), the matrix 
equation may be inverted directly. However, the kernel G i} is 
ill-conditioned and sensitive both to error and numerical 
rounding. The art therefore, is to provide methods of inver- 
sion which can deal with noise. 

An obvious approach is the method of least squares. 
Perhaps the most popular such technique is that due to Phil- 
lips 1465] and Tworoey (466]. However, it has the disadvan- 
tage that the operator has to adjust an unknown Lagrangian 
multiplier to achieve an appropriately smooth output curve. 
Much is left to personal judgement, although Amato et al. 
(467] discuss ways of selecting the smoothing parameter. In 
a further development, Shifrin and Zolotov [4681 have 
suggested a means of choosing the most appropriate regu- 
larization parameter. 

A modified form of the Phillips-Twomey method was 
used by Ramachandran and Leith (469). Their results were 
satisfactory for 0.34 £ A =s 0.94 am and 0.3 £ D ^ 
2.5 /im. Ramachandran et al. [470] have further developed 
the method by the use of computational tomography to 
enable concentration and size distribution to be obtained 
as a function of position in an aerosol cloud. 

Other methods of direct inversion have been discussed by 
Bertero et al. [471, 472], Ben-David et al. [473, 474]. Viera 
and Box [475-477], Curry [4781, Shimizu and Ishimaru 
[479], Box and Viera [480], Wolfenberger and Seinfeld 
[481], Brianzi and Fontini [482], Schnablegger and Glatter 
[483], Box et al. 1484], Dellago and Horvath [485. 486], 
Ferri et al. [487], Dubovik et al. [488], Gail [489], Bryant 
and Thomas [490] and Ramachandran and Kandiikar [491]. 
An interesting variant due to Ligon et al. [492] employs a 
Monte Carlo technique in which sizes are chosen at random 
and the sum of the squares of the difference between the 
predicted and measured scattering parameters is minimized. 
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Benayahu et al. [493] have attempted direct inversion in the 
presence of double scattering events. The influence of 
refractive index on the inversion problem in the anomalous 
diffraction regime has been investigated by Oshchcpkov and 
Dubovik [494]. 

Recently there seems to have been some renewed interest 
in fitting the experimental data using presupposed trial func- 
tions for the size distribution. Sun and Lewis [495] assumed 
a log-normal distribution and adjusted the parameters to 
minimize the sum of the squares using a simplex deconvo- 
lution. Amato et al. [496] compared the Twomey method 
with a multi-modal log-normal distribution and found 
reasonable agreement. 

A lengthy review has been given by Bates et al. [497]. 
They propose an approach which saves computer time by 
sacrificing accuracy. Nonetheless, the accuracy is reason- 
able. 



10. Measurement of refractive index 

In order to fully characterize particles, their refractive 
index should be measured as well as size and shape. Not 
only would this enable panicle types to be identified, but 
proposals are being made to measure temperature directly 
through the variation of refractive index. 

Addition of refractive index increases the number of vari- 
ables to be measured. For spheres this becomes three — size 
and the real and imaginary parts of the refractive index. Any 
one scattering measurement is a function of all three, 
although close to the forward direction size is by far the 
dominant parameter. Thus a single measurement can predict 
a range of possible combinations of the variables lying on 
one or more surfaces in the variable space. Multiple 
measurements, in principle, can pinpoint the true values 
by noting where all the surfaces intersect. This approach 
was attempted by Avery and Jones [415] for coal, though 
this was unsuccessful because of the particle irregularity. 
Similar problems were encountered by Bhattacharya et al. 
[498]. 

For more nearly spherical particles, intersecting contours 
obtained by the visibility method were successful for non- 
absorbing particles [499]. Later measurement of the full 
complex refractive index of absorbing liquid drops was 
demonstrated by Card and Jones [500], again from the visi- 
bility method, although a total of four scattering angles was 
needed. Takeda et al. [501 ] have also used contour matching 
and employed two wavelengths and two angles to provide 
four measurements. However, only the real part was 
obtained and dispersion was neglected. Methods involving 
minimizing the difference between measured and scattered 
light parameters have also been proposed by Chae and Lee 
[502] and Schnablegger and Glatter [503], 

The possibility of material recognition via refractive 
index measurement has been proposed by Naqwi et al. 
(423, 504]. The suggested method involves measuring 
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fig. 33. Magnitude squared of the Fourier transform of the angular 
scattering diagram near the rainbow. (After van Beeck and Ricth- 
muller (514].) 

both phase and frequency of the scattered light. Extra 
detectors are required and the refractive index comes from 
the ratios of the phases or frequencies. 

Another interesting possibility for the measurement of 
refractive index and size for the case in which the particle 
diameter is much larger than the beam width, is the so-called 
dual burst technique [505]. The principle behind this can be 
understood by considering the behaviour of a single ray as it 
is traversed by a drop, as in Fig. 32. At any one angle there 
are two pulses* one due to reflection off the outer surface and 
one due to refraction. For a given panicle size the time delay 
is a known function of the refractive index. Also, if there is 
absorption, the height of the pulse due to refraction is 
reduced. Thus it is possible in principle to measure the 
full complex refractive index. The authors have performed 
a full analysis for the case of a narrow Gaussian beam, and 
have also investigated the prospect of determining inhomo- 
geneity of refractive index by examining concentric spheres 
1506]. 

A specific suggestion for refractive index measurement 
with a view to remote temperature sensing has come from 
Maxsoli et a). [507]. They propose that the ratio of scattered 
intensities at 33° and 60* be measured by using horizontal 
polarization. This displays a strong monotonic dependence 
on refractive index, although there is some sensitivity to 
size. The technique only applies to homogeneous spheres. 

Clapp et al. [508] measured the absorption spectrum of an 
aerosol of water ice in order to obtain the imaginary part of 
the refractive index. Hie real part was found by use of the 
Kramers-Kronig [191] relationship. They quote the full 
complex refractive index of ice as a function of both 
temperature and wavelength over the range 2.5 pjn to 
12.5 tun. 

A different approach has been discussed by Taubenblatt 
and Batcheider |509], essentially employing interferometry. 
The scattered light is mixed with a reference wave so that 
both amplitude and phase can be measured. A plot of the 
extinction against phase shift is found to be sensitive to 



refractive index. Even for non-spherical particles it is 
claimed that there is sufficient information for particle 
discrimination. 

An obvious candidate for the measurement of the refrac- 
tive index of large transparent spheres is to determine the 
angular position of the rainbow. For spheres which are suffi- 
ciently large, this is independent of size. Massoli et al. [507] 
showed that at optical wavelengths the diameter needed to 
be greater than about 60 urn 

The position of the primary rainbow is found by geo- 
metrical optics as the angle^of minimum deviation of a ray 
which has undergone two interna) reflections. There are 
higher order rainbows associated with multiple internal 
reflections. The positions of the rainbows are given by 
Bohren and Huffman [4] as 



cos e m 



- 



m l - 1 
n(n + 2)' 



ft — w 1 y 2* 3, I 



where m is the refractive index. 

An attempt to apply the rainbow method to the measure- 
ment of refractive index of burning fuel drops, and through 
this to obtain the temperature, was made by Roth et al. 
[510]. However, it was soon pointed out thai such a drop 
would have a temperature gradient within it and that this 
may substantially influence the result of a rainbow measure- 
ment. Schneider et al. [51 1] performed an analysis on this 
and concluded that errors in the measured temperature 
would be significant. Kai et al. (25] used the stratified sphere 
model to indicate the ways in which this position fluctuates 
and Anders et al. [512] used geometrical optics to the same 
end. The latter authors suggest that if the surface tempera- 
ture could be measured independently, then rainbows could 
be used to find the internal temperature gradient. 

Inhomogeneiry has also been investigated by Lock et al. 
[513]. In this case they examined the behaviour of the rain- 
bow for a coated sphere as the thickness of the coating 
varied. For a coating thickness 6 and a drop radius A, the 
intensity of the rainbow oscillates with thickness provided 
67* I0~\ For S/R * !0" 2 the primary rainbow breaks 
into two separate rainbows. 

An interesting technique to measure particle size and to 
indicate shape has been described by van Beeck and Rteth- 
muller [5141. They take the Fourier transform of the scatter- 
ing pattern for angles equal to and greater than the rainbow 
angle. The amplitude spectrum produced has a series of 
peaks, as indicated in Fig. 33. The positions of the various 
peaks are given by 



sinr / 

2.3794 { A 2 cos 
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where 



W - 1 

The size can be determined cither from F, or F 2 . If the two 
sizes are not equal, the panicle is. not spherical. 

Some interesting phenomena occur when the beamwidth 
is of the order or less than the diameter of the particles. For 
example, Lock (531 demonstrated that a narrow Gaussian 
beam incident near the edge of a particle can produce rain- 
bows up to ninth order. Chan and Lee [515] observed up to 
32 rainbows. They showed that a plot of angular position 
against order could be used to obtain refractive index. 

11. Tomography 

With the advent of cheap computing power a new method 
has become available for the reconstruction of the structure 
of objects. This is tomography, in which the object is 
subjected to a test and the response is observed remotely. 
The test may be the application of a magnetic field, a flux of 
neutrons or X-rays, as in medical tomography. In scattering 
a cloud of particles is illuminated by electromagnetic radia- 
tion and the scattered intensity is found. The iUurninating 
field may be incident on all points of the object simulta- 
neously, or it may be scanned across. Similarly, the remote 
response may be measured at all points simultaneously* or 
the distribution may be built up over time. 

The relationship between the detected pattern and the 
structure of the object leads to a complicated set of mathe- 
matical equations. It is for this reason that significant 
computing power is required. 

An early example of the application of tomography to a 
scattering system was the application to a sooting flame by 
Jagoda et aL [516]. The problem was simplified by assuming 
axial symmetry. In this method the cloud is divided into a 
number of concentric rings each sufficiently narrow that the 
concentration within it may be considered uniform, as in 
Fig. 34. The laser beam first passes only through the outer 
ring along path 1. Since the properties are uniform, the size 
and concentration are easily determined The beam next 
follows path 2. It suffers known extinction as it crosses 
ring I on the way in and out. The extinction due to ring 2 
alone is then easily found. Similarly, the light scattered out 




Fig. 34. Illustration of the subdivision of an axially symmetric 
system into a set of concentric circles. 



of ring 2 suffers known extinction as it emerges through ring 
1. Thus the size and concentration in ring 2 are found. The 
process is repeated until the central zone is reached. 

A more genera) form of this technique replaces the scan- 
ning of the beam and single detector by an expanded beam 
illuminating the full diameter of the flame and an array of 
detectors. There is one detector for each ring, as above. For 
systems which are not axi symmetric the beam -detector 
assembly can be rotated around the flame and measurements 
made at a number of angles, as in medical tomography. 

More recent applications to scattering include work by 
Sivanathu and Gore [517] and Menguc; and Dutta (518]. In 
both cases axial symmetry was assumed. In the latter, two 
beams were used uniformly spaced about the centre of a 
sooting flame, the comparison of the two yielding an absorp- 
tion coeflicient profile. Bates et al. [519] applied Fourier 
transform Hadaraard tomography to sooting flames and 
Drescher et al. [520] measured gas concentrations in air. 
Tomography has been applied to aerosols by Ramachandran 
et al. [470] and to sprays by Oberle and Ashgriz [521] using 
back scatter. A Doppler tomography method has been 
devised by Wang et al. [522] to map fluid flow velocity. 



12. Conclusion 

It is now of the order of 100 years since Lorenz ( 16). Mie 
[17], Rayleigh [12] and Tyndall [523] laid the foundations 
of light scattering. Since then the technique has developed 
into one of the most powerful diagnostic tools for probing 
the properties of particulate matter. 

The primary focus of attention for panicle characteriza- 
tion has been the measurement of size. Traditionally this 
was based on the assumption of a sphere, because of 
theoretical restrictions, and considerable effort went into 
developing techniques which were insensitive to shape 
and (often unknown) refractive index. These were largely 
successful, and more recent developments have been 
concerned with extending the dynamic range of the instru- 
ments, with regard to both the size and the range of particle 
concentrations over which they may be employed. The 
measurement of size is now largely a matter of course, 
although improvements in instrumentation continue. 

Interest has now begun to switch to the determination of 
other characteristics, such as shape, structure, composition 
and the distribution of material within a particle cloud. From 
small beginnings about a decade ago, the rate of publica- 
tions in this area is accelerating rapidly. Some of the areas of 
particular interest to combustion include the determination 
of drop temperature in a spray, the structure of soot agglom- 
erates, the measurement of residual carbon in Ay ash and 
tomography to reveal the structure of a burning spray. 

On the theoretical side, there have been massive improve- 
ments in numerical methods for the calculation of the scat- 
tering properties of non-spherical and non-uniform 
particles, and in inversion techniques for the direct determi- 
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nation of particle properties from light scattering data. The 
new mathematics of fractal analysis has had considerable 
impact on the understanding of agglomerates of small 
particles. 

Even after a century light scattering continues to be a very 
vigorous Meld of study. This is evident from the large and 
growing number of publications appearing each year. There 
is every reason to believe that the field will continue to grow 
with many exciting developments to come. The growing 
power and speed of computers will inevitably enable the 
full potential of numerical methods to be realized, so that 
the properties of particles of complex shape and structure 
can be calculated for practical values of size and retractive 
index. They will also assist in the study of the statistical 
nature of random shapes and multiple scattering interac- 
tions, and facilitate direct inversion procedures and tomo- 
graphy. No doubt the extra information generated will lead 
to new measurement techniques for a wide range of particle 
properties, especially in those areas where interest has 
emerged relatively recently such as irregularity, refractive 
index, composition and structure. General technological 
development such as new light sources, improved optical 
fibres, solid-state devices and faster electronics will also 
have an impact, leading to enhanced and novel instrumenta- 
tion covering a wider range of operational situations. In the 
combustion context, this would include the ability to probe 
hostile environments with increased ease and confidence. 
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